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ABSTRACT 


The mathematical model of geopolitics is a conditional name for 
several models, which are naturally connected and act as 
accompaniment to the main theme — geopolitics. All constructed 
work models made the transition to computing experiment, the 
results of which are given and discussed. The central concept of the 
mathematical model of geopolitics is introduced — the capacity of 
the habitat. Geopolitics is the climate, the relief, the logistics features 
of global commodity flows, the geopolitical confrontation in terms of 
the "sea-continent", i.e. all that constitutes the material complex of 
living conditions of the inhabitants of the Earth. This complex, to a 
large extent, mediates the population's behavior from the political 
point of view. The author does not adhere to the position of natural 
determinism, which acts in the form of geopolitics, but tries to 
delineate the scope of the manifestation of geopolitics in real politics. 
The paper provides clarification and generalization of the generally 
accepted geopolitical classification of territories in terms of 
orientation and positioning either at sea or on the continent. The 
mathematical model of transport costs minimization for an arbitraiy 
number of points acting in the form of logistic centers is formulated. 

Keywords: the capacity of the habitat, the average annual temperature 
and moisture, relief, correlation analysis, cutting territories, the 
geopolitical classification of points [territories], traffic, nonlinear 
optimization. 
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INTRODUCTION 


In a number of works by the author [1, 2] the "Normative model of global history" was 
formulated. At the heart of this model lies the representation of global political integrity as 
a set of atomic units, called geopatom (GEOPolitical ATOMs). It seems that the current 
political structure of the territories of human habitation on Earth, to a large extent, can be 
explained by combining a specially prepared smaller territories — geopatoms. In 
connection with this hypothesis there is the task of building a set of geopatoms, covering 
the entire territory of the Earth. 

For adequate cutting of the Earth's territory into individual geopatoms, it is necessary to 
take into account the climatic features of the territories, as well as the relief. In addition, it 
is also necessary to take into account global traffic, which significantly changes the 
configuration of the location and shape of the geopatoms. 

To take into account for climatic peculiarities of living in these territories the concept of 
"capacity of the habitat" will be formulated. Will be built a suitable density function of the 
capacity of the habitat, which will include the whole complex of anthropomorphic 
requirements to the environment in terms of climate. The areas of residence of people will 
be classified in terms of "[not) highly - (not] favorable", i.e. taking into account height 
above sea level and climatic features. 

In terms of global traffic modelling, the central geopolitical confrontation between the "sea" 
and the "continent" will be particularly acute, as warned by the classics of geopolitics. The 
formulation of the algorithm and the calculation of global traffic flows with its help will 
clarify and summarize the classification of all territories in the conventional geopolitical 
categories. 
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SECTION I. 

THE CAPACITY OF THE HABITAT 


We will construct a mathematical model of what is called differently in different 
disciplines: the capacity of the habitat, the living space etc. We will be interested in the 
construction of the measure of the volume the capacity of the habitat by taking into account 
various kinds of the most important natural components that accompany people's lives. 
Such components include, primarily, the annual averages of temperature and precipitation 
on the Earth's surface. 


A spatial scan in the model bind with latitude cp and longitude X of the Earth's surface, 
considering that [-;r/2, .7r/2], X&{-7r,7^ — in radians or ^ 3 e[—90°, 90°], 

/le[—180°, 180°] — in degrees. To select a suitable formula for the connection of the 
density of the capacity of the habitat, p with the average annual temperature Tannuai 
the precipitation Qannuai we draw real data available for the entire surface of the Earth. 


We introduce comfortable temperatures from the point of view of people's living, Tcomfort 
and precipitation, Qcomfort according to the weighted average procedure, choosing the 
density of people's living on the surface of the Earth as a weight, Pdensity this case, we 
can write down the following general formula for calculating Tcomfort Qcomfort- 


rp _ JJ ^annualPdensity (pdcpdA _ JJ QannualPdensity ri'\ 

comfort JJ Pdensity cos ’ comfort JJ pdensity cos (pd<pdA ^ 


where ^annual '^annual ^annual QannualC^P'Pdensity Pdensity 

corresponding latitude and longitude functions. The integration in (1) is assumed 
throughout the globe without taking into account the water surface. 


Let's choose a spatial resolution of 0,5° degrees and consider the specific data for the 
average annual temperature, precipitation and population density. In this case, the data 
will appear in the form of matrices of size 360x720, where 360 is the number of nodes of a 
uniform grid in latitude and 720 — longitude. Grids in latitude and longitude in radians can 
be represented as cpi = ^[“89,75 + 0,5(i — l)],i = 1, ...,360; Xj = ^[—179,75 + 
0,50 -1)],; = 1,2, ...,720. 

Load the data on annual average temperature and precipitation from the site [3]. We 
choose the version of archive 2.02, which is dated July 2001. A description of the archive 
data is available on the website [4]. Data on the density of the population of the Earth we 
take from the site NASA Earth Observations [5]. 


Replace the integration in (1) by summation over our grid, then 


rp EiJ rannual,i,jPdensity,i,j COS (pj ^ Eij 0annual,i,jPdensity,i,j cos (pj 

^comfort = v~Tr ^ ”77777 ' Vcomfort = y7~77 ^ r”77777 ’ 

Li,] Pdensity,1,7 Pdensity,1,7 


( 2 ) 
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where the summation in ( 2 ) is carried out over those grid nodes that do not belong to the 
water surface. 

According to the formula (2) the calculation was carried out with data from the archive 
2.02 and NASA data on population density, it was found that Tcomfort = 18,5775°C; 
^comfort 1200,2 mm. 

The weighted average standard deviations from comfort values were also calculated, they 
turned out to be equal to Tcomfortstd = 7,9084°C and Qcomfortstd = 722,5299 mm 
respectively for temperature and precipitation. Such noticeable standard deviations 
indicate that the proposed dependence of the capacitance of the medium p from deviations 
from the comfort values in the form of power expressions such | T/rcomfort ~ 1 |“ 

I Q/Qcomfort “ 1 |^ should be small compared with 1 values of degrees or and p. 

After entering comfort values of temperature and precipitation, we choose for further 
analysis the following formula for describing the density of capacity of the habitat: 


where a-^, 02 , a^, a,p — while undefined constant. 

In the representation (3) the entire set of requirements for the function p is focused. The 
function describing the capacity of the medium^, p must be nonnegative, depend on 
deviations from comfort values that are easily modified if an additional parameter 
describing the capacity of the medium is taken into account. 

The uncertain parameters a-i, a 2 , entering in (3) are found by minimizing the sum of 
squares of deviations of the inverse function of the medium's capacity, 1 /p from the 
inverse population density function, 1 /Pdensity It turned out, that precisely in this case the 
correlation between the inverse of the capacity of the medium and the inverse of the 
population density, Sdensity = 1/Pdensity is maximized. The value of the Sdensity a clear 
geographical sense — it is the number of square kilometers per person in the region, if the 
Pdensity is measured in the number of people per square kilometer. 


Taking into account (1), we write down a functional D whose minimum provides a search 
for the parameters a^, a 2 , for fixed values of a, P: 


D — Z) U 2 , U 3 ) 


flip H<P,A)-Pd^nsity]^cos(pd^dA 
JJ cos (pdcpdA. 


min D, (4) 


where integration is assumed over the entire land, without taking into account the water 
surface. 


1 Means the same as "the capacity of the habitat” or "the capacity of environment”. 
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Taking (2) into account, we rewrite the functional (4) in finite-difference form, proceeding, 
as above, from the resolution of 0,5°. In this case, up to the numerical coefficient, we obtain 

^ D(^Cl±, 0,2) ^ 3 } — jj^i,j\_PiJ Pdensity,i,y] COS (pi, (5) 

where i = 1, ...,360; j = 1, ...,720; N — the number of grid nodes on the land surface, and 
the summation extends only to those grid nodes that fall on the land. 

To find the minimum of the functional in the form of (5) the problem of gradient descent 
was solved 


dttfc 30(01,02,03) , 1 o o 

— = -^-,/C = 1,2,3, (6) 

dt duk 

where t is some conditional argument, herewith 

do 2 -1 -1 -1 

^ ^Pi.i -Pdensity,u)cos<Pi, 

2 r? —1 —1 —1 

~ J^^l,j\Qi,j/Qcomfort~^V PiJ (PiJ ~Pdensity,i,;) ‘Pi- 

The system of differential equations (6) can be solved using one of the standard solvers, for 
example, the MATLAB package. The solution was carried out until the values a^, a 2 , were 
established, i.e. when it is possible to clearly identify the presence of the limits of the values 
tti, a 2 , for t -l-oo. 

Table 1. List of values of the parameters included in the formula for the capacity of 

the habitat 


# 

^comfort 

Ocomfort 


«2 

«3 

a 


Corr 

1 

18,5775 

1200,2 

0,7156 

-0,0649 

-0,0426 

2 

2 

0,5758 

2 

18,5775 

1200,2 

1,0227 

-0,2767 

-0,2991 

1 

1 

0,6329 

3 

18,5775 

1200,2 

1,6592 

-0,6627 

-0,6656 

0,5 

0,5 

0,6581 

4 

18,5775 

1200,2 

6,2794 

-3,2566 

-2,7573 

0,1 

0,1 

0,6590 

5 

18,5775 

1200,2 

3,4189 

-1,6624 

-1,4789 

0,2 

0,2 

0,6611 


Table 1 contains a list of the values of the parameters entering into the formula (3) for the 
medium capacity, obtained taking into account the solution of problems (2), (5), (6). In 
addition. Table 1 shows the corresponding values of the correlation coefficient, Corr, which 
establishes the relationship between the sets [plJ} and {Pdensity.i,;}' ^ = 1, ...,360; j = 
1, ...,720. Taking into account the features of the geometry on the sphere, the correlation 
was calculated by the formula: 


Corr = 


P ^^’-Pdensity.y Pdensity^“Pi 


P ^i'JN^iJ^^density,i,J ^density^^ 


( 7 ) 
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where 


P PiJ COS (Pi, PfJensity Pdensity,i,j cos (pi- 

Figure 1 shows typical samples of the dependence of the modulus of the gradient G = 

I 2 2 2 

J(^) "*■ (^) ■*" (^) value of the functional (5), D on the argument t. The 

graphs in Figure 1 demonstrate the obvious convergence of the modulus of the gradient G 
to zero (the left graph), and the value of the functional D to some limiting, minimum value 
(right graph). This means that the set of parameters a^,a 2 ,a^ converges to its limiting 
values when + go. The other parameter values are chosen according to point #5 of the 
Table 1, when there was a maximum correlation value. 


Dynamics of the dependence of the gradient 



The dynamics of the dependence of the value 



Figure 1. D 5 mamics of modulus gradient G (left graph) and the functional values D (right graph) on 

the argument t 


The negative values of the signs of the parameters 02 , in Table 1 meet our expectations. 
The values of parameters a and (3 were subjected to variation. It was found that the 
maximum correlation value, excluding a small variation is in the neighborhood of values of 
degrees a = 0,2 and = 0,2. Note that if in the formula for the capacity of the medium (3) 
the precipitation is not taken into account (a 3 = 0), the correlation coefficient at the 
parameter value from point #5 of Table 1 will be 0,6341. The resulting value is slightly less 
than the number 0,6611, i.e. precipitation component provides its contribution. Finally, in 
the absence of temperature in formula (3) { 0.2 = 0) the correlation value was 0,5587, 
which is noticeably less than 0,6611. These observations mean that the combined use of 
temperature and precipitation indicators increases the correlation between the values of 
the inverse capacity function of the habitat and the inverse function of population density. 
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Thus, it can be concluded that approximately two-thirds of the correlation between the 
inverse function of the habitat capacity calculated by formula (3) and the inverse 
population density function falls on the habitat in terms of temperature and precipitation. 

Let us compare visually the maps of the inverse function of the density of the capacity of 
the habitat and the inverse function of the population density. The required maps will be 
drawn using level lines. Figure 2 shows the result. On the left map of Figure 2 shows the 
level line of the inverse function of the density of the medium's capacitance, calculated by 
the formula (3) with parameters from point #5 of the Table 1. In total, 12 level lines with 
equidistant values are constructed, which show the range of variation of the invers capacity 
[0,0054;1,1268]. On the right the map Figure 2 shows the inverse function of population 
density. Also constructed are 12 level lines with equidistant values, which represent the 
range of variation in the population inverse density [0,0001;0,9664]. 

In Figure 2 more blue up to violet flowers with small values of the inverse function of the 
density of the habitat capacity and inverse population density function are marked. There 
is a clear correlation between the locations of the blue and violet spots on the left and right 
maps of Figure 2. Note that these places meet the most comfortable place for human living 
conditions. There are significant differences on all continents. To the right of each of the 
maps is a palette that establishes a connection with the numerical value of the 
corresponding indicator. In addition, on each of the maps has a coastline in the form of a 
red line. 
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The level lines of the inverse function of 
population density 


OCA 


0 100 200 300 400 500 600 700 


The level lines of the inverse density function 
of the habitat capacity 




0 100 200 300 400 500 600 700 


Figure 2. Maps of level lines of the inverse density function of the habitat capacity 
(left) and the inverse population density function (on the right) 

Let's move from the inverse density values of the habitat capacity and density of population 
to the direct values. Compare the linear Pearson correlation coefficient and rank 
correlation Kendall and Spearman correlation for the density capacity of the habitat and 
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density of population. To do this, we transform the matrices Pij^cos (pi and 
P density,ijyj cos (pi into vectors and find the corresponding correlation coefficients between 
these vectors. Table 2 summarizes the results. 


Table 2. The correlation coefficients of Pearson, Kendall, and Spearman 
between the density of habitat capacity and the population density 


Correlation coefficient 

Pearson 

Kendall 

Spearman 

The value of the correlation coefficient 

0,1740 

0,7582 

0,9034 


Table 2 shows how large are the rank correlation coefficients compared with the 
conventional linear correlation coefficient. For this reason, in the formula (7) the linear 
correlation between the inverse density function of the capacity of the habitat and density 
of population was considered. For the inverse values, the linear correlation coefficient was 
found to be 0,6611, which is much larger than 0,1740. 

As a result of the above constructions and calculations, we choose the function (3) for the 
function of the capacitance of the habitat of people, which reflects the relationship with 
such characteristics as temperature and precipitation, or more precisely, with deviations in 
temperature and precipitation from the corresponding comfort values. 
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SECTION II. 

THE DISTRIBUTION CAPACITY OF HABITAT IN THE STATES 


Taking into account the functional representation for the habitat capacity (3), we will write 
the formula for calculating the global capacity of the habitat: 

U j^^dAcos(p ■ p(T((p,A),Q((p,A)). (8) 

In the beginning, we will find the center of concentration or the center of the scattering of 
the capacity of the habitat and compare it with the center of the dispersion of the 
population of the Earth. Let us turn to a discrete representation of all integral indicators, 
starting from a grid in latitude and longitude, the size of 360x720. We calculate the center 
of population dispersion according to formulas: 

_ SiJ <PiPdensity,i,y cos _ 2i,y-^y'Pdensity.iJ cos (p, 

<Ppopulation - 2yPdensityX;COS<Pi ' '^PoPulation - PdensityX; cos 

where the summation extends to those pairs (^j,Ay), i = l, ...,360; y = 1, ...,720, that 
belong to the land. The population density of Pdencity be the same as in the previous 
section. 

Similarly to formulas (9), we define the expressions for calculating the scattering center of 
the capacity of the habitat: 


_'Zl,j(PiPij cos <PI cos (Pi 

^Pcapacity ~ v ' ^capacity ~ v ^ 

^ hi,j Pi,j cos (pi c i,i j Pi j cos <Pi 

where pij, i = 1, ...,360; j = 1, ...,720 is the density of the habitat capacity (3) with the 
values of the parameters from point #5 of Table 1. 

Figure 3,a shows the total count of population dispersion centers, (^populations-^population) 
and habitat capacity, (^Pcapacity< ^capacity)- b turned out that the center of population 
dispersion has coordinates (23,78 N; 51,85 E) — Arabian Peninsula, Saudi Arabia, salines. 
The center for the dispersion of habitat capacity has coordinates (13,59 N; 19,03 E) — 
Central Africa, Central Chad, east of Ndjamena. 

Note that in the model, both the dimension and the absolute value of the habitat capacity 
are not significant. Only the spatial sweep is important, i.e. the relative values of the 
capacity of the habitat or individual regions, or the density of the medium of individual 
points. 

As an example of calculating the global capacity of the environment (8) let us consider the 
world political map with borders of individual States and territories and calculate the 
capacity of the habitat for each of the States and territories. According to the UN data for 
2008, the borders of 246 States and territories, data on which are presented, for example. 
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on the website [6], are marked. Figure 3,b shows a graphic illustration of all the boundaries 
of states and territories available for 2008. 


The centers of scattering of population and habitat capacit> 



Figure 3,a. Population scattering centers 
(Arabian Peninsula), and the capacity of 
habitat (Central Africa) 


The borders of 246 states and territories of the w orld 



Figure 3,b. The borders of 246 states 
and territories of the world 


Define the region occupied by the k-th State or territory, with the symbol In this case, 
the global capacity of the habitat (8), taking into account the selected marking of the 
Earth's surface, can be rewritten as: 

U = 'Zlt\uk,Uk = ff ^^d(pdA cos (p ■ p(T,Q). (11) 


To calculate the capacity of the habitat of all States and territories Ui^, k = 1, ...,246 we 
replace the double integrals in (11) with approximate sums, starting from a resolution of 
half a degree in latitude and longitude. In this case, the second equation in (11) can be 
rewritten in the form: 


2 

Uk = (jk) IcpajSWk^OSCPiPij, ( 12 ) 

where i = 1, ...360, j = 1, ...,720, and the summation in (12) extends to those points 
{(Pi, Xj), that enter the region 

A computational experiment was carried out to calculate the set of capacities of the habitat 
of States and territories (12), while the matrix pi j was calculated from formula (3) with 
parameters from point #5 of Table 1. 

Figure 4 shows the final graph, which marked some countries and territories, deployed in 
terms of habitat capacity. In particular, the top five in terms of capacity of the habitat in 
descending order is as follows: Russia, the United States, Brazil, China, Australia with the 
values of 0,6027; 0,5818; 0,5261; 0,5061; 0,4968 respectively. 

The global capacity of the habitat U was 7,6438, with an average value of — 0,0311, and a 
standard deviation of — 0,0863. 
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The distribution of the States and territories in terms of the habitat capacity 
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Figure 4. The distribution of the States and territories in terms of the habitat 

capacity 


Let us study the distribution of the capacity of the habitat per person. We introduce a 
suitable variable in the following form: 


ffc ~ tt;c/Ppopulation,fc< (13) 

where Ppopuiation./c — population of the k-th State or territory, k = 1 ,... ,246. 

The distribution of the States and territories in terms of the specific habitat capacity 
1 ^Greenland 
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Figure 5. The distribution of the States and territories in terms of the specific habitat 

capacity 
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The variable introduced in (13), denotes the volume of the habitat per person within 
the k-th State or territory or the specific value of the volume of the habitat. As specific 
values of the population were selected data for 2005, which are given in the same data 
source as the data of the borders of 246 States and territories. 

The performed calculation of the set of indicators k = 1 ,... ,246 showed that the winner 
is Greenland. Figure 5 shows a graph similar to that in Figure 4, where some values of the 
set of specific volumes of the medium k = 1, ...,246, normalized to the maximum 
specific volume Vcreeniand' plotted along the ordinate. 


Tables. The ratio of the specific volume of the Russian habitat to 
_ the specific volumes of a number of other countries _ 


Russia/Finland 

Russia/United Kingdom 

Russia/France 

Russia/Germany 

1,5416 

17,6738 

7,5665 

16,6839 

Russia/Poland 

Russia/Ukraine 

Russia/Turkey 

Russia/Kazakhstan 

10,4924 

6,8427 

6,7931 

0,5802 

Russia/China 

Russia/Japan 

Russia/USA 


10,8619 

23,3693 

2,1576 



Consider the ratio of the specific capacity of Russia's habitat in relation to a number of 
other States. Table 3 shows some characteristic list of relations. 


Arrows of abnormal values of the ratio of the specific habitat capacity of 
the Russian Federation and specific capacities of other countries 



Figure 6. Arrows of abnormal values of the ratio of the specific habitat capacity of the 
Russian Federation and specific capacities of other countries 

According to Table 3, the largest values of the ratio of the specific volume of the Russia 
habitat to the specific volumes of other countries fall on the three countries traditionally 
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included in the historical orbit of the Russian Federation: The United Kingdom, Germany 
and Japan. 

We will construct a map on which the arrows indicate the transitions from those countries 
which, on the one hand, have a sufficiently high potential for the capacity of the habitat 
(Ufc > 0,02 max Uj), on the other hand, the ratio of the specific habitat capacity of the 

l<i<246 

Russian Federation to the specific capacity of another country is greater than a certain 
threshold value (^Russia/t^fc > ^'5). Thus, we take into account only those countries whose 
habitat capacity exceeds 2% of the capacity of the largest state, i.e. the Russian Federation. 
As the threshold value of the ratio of the specific habitat capacity of the Russian Federation 
and the specific capacities of other countries, we choose the value 7,5, that is, a sufficiently 
high value. 

In Figure 6 shows a map with the desired arrows. With the selected thresholds for the 
capacities of the habitat and the specific capacities of individual countries, there were 13 
arrows that oppose 13 countries with the Russian Federation. Part of the countries shown 
in Figure 6 also appear in Table 3. Note that all 13 arrows are directed to the Russian 
Federation. 
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SECTION III. 

CAPACITY OF HABITAT (POPULATION DENSITY) AND RELIEF 


Using the relevant data, we will study the density of the habitat capacity (population 
density) in connection with the relief of the Earth's surface^. The relief data is represented 
in latitude, longitude coordinates with a resolution of 0,5° in the form of a matrix Z^j, i = 
1,...,360;; = 1,...,720. 

We define a set of average heights in different senses: the average height of the land, Z; 
weighted average height, Zpopujation^ where the population density is chosen as the weight; 
as well as the average weighted population density of 50% of the population, 
^population, 50 %. where Z < Z* and Pdensity > P*> where Z*, p* — are some nonnegative 
constants. The averages listed above are found according to the following formulas: 

y _ 2i,jCOS <P( - _ Ei,j 2i,yPdensity,i,y COS <P( 

~ hjcoscp, ’ ^population - Pdensity,y cos.p, ' 

The formula for calculating ^popuiation, 50 % is similar to the second formula in (14), provided 
that the summation extends over those pairs of indices (i,)), for which the following 
inequalities are true: Zij < Z* and Pdensity,!,; > V*- The threshold values Z*, p* were selected 
from those considerations that the proportion of the population living under the conditions 
Z < Z» and Pdensity > V* was equal to 50%. The result of the calculations of the average 
values and parameters is given in Table 4. 


Table 4. A set of average heights and suitable parameter values 


z 

^population 

^population, 50% 

z. 

P* 

804,9 m 

442,9 m 

61,3 m 

187,5 m 

2 people/sq.km 


According to Table 4, the average height of land is 804,9 meters. Weighted average height 
of the land, taking into account the population density, roughly speaking, is half the size, 
more precisely, 442,9 m. Finally, 50% of the population lives within the height of 187,5 m, 
while the average weighted height of the settlement, taking into account the density of the 
population, is 61,3 m. The last mean value characterizes the well-known fact that half of the 
population is concentrated in the vicinity of the coastline of the world ocean. 

Figure 7 shows two maps of settlement of 50% of the population concentrated near the 
coastal edge and provided that the density of the population exceeds the selected threshold 
p». Figure 7,a shows the coastline (red) and region of settlement population (toned in 
blue). On the map Figure 7,b with different colors in a separate closed polygons that meet 
the conditions that people live at altitudes not exceeding the value of Z„ with a density of 
population exceeding the value of p„ are shown. In this calculation, the total number of such 


2 The relief data was taken from the Mapping Toolbox of the MATLAB package. 
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polygons was 932, while the shoreline is not marked, but it is readily apparent in the 
position of the polygons. 



Figure 7,a. Map of the resettlement of 
50% of the Earth's population (overall 
perspective) 


Map of the resettlement of 50% of the Earth’s population 





Figure 7,a. Map of the resettlement of 
50% of the Earth's population (multi- 
polygonal perspective) 


The above considerations on the relationship between the relief and the population density 
of the Earth's population speak of their complex interrelationship. In particular, it is not 
possible to introduce a comfortable value of the altitude, like comfortable temperature and 
precipitation. 


We will solve the problem of identifying areas on which 50% of the total capacity of the 
habitat U is concentrated. These territories can be found from the condition that the 
heights do not exceed a certain threshold value Zp, and the density of the medium capacity 
exceeds a certain threshold value Pp. As the density of the capacity of the habitat, we take 
function (3) with parameters from point #5 of Table 1. We define the weighted average 
height Zcapacity< where the density of the medium capacity p, appears as the weight, and 
also the average weighted average over the territories, Zcapacity, 5 o%< where 50% of the total 
capacity of the medium is concentrated. Similarly to (14) we write the formula for 
calculating the weighted average height: 

7 _ 'Ll,] ZijPij cos (pI 

capacity Y.i.jPi,, cos 

where the summation in (15) propagates over those pairs (i,;), for which Zjj > 0. 


Table 5. A set of weighted average heights and suitable thresholds for 


parameters 


aking into account the density of the mediu m's capacity 


7 

^ capacity 

^capacity, 50 


Vp 

U 

714,5 m 

312,1m 

820 m 

1,82 

7,7039 


Table 5 shows a set of parameter values, which is the solution for the problem formulated 
above. 
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The parameters given in Table 5 allow us to find the areas where 50% of the habitat 
capacity is concentrated. These are the territories that are within the height Zp = 820 m, 
and the values of the density of the habitat capacity exceed the value Pp = 1,82. In Figure 8 
maps are constructed that illustrate the positioning of the 741 territories that fall within 
these limits. On the left map of Figure 8 the search area is shaded one color and on the right 
the map different colors. 

Map of territories that carry 50% of Map of territories that carry 50% of 




Figure 8. Map of territories that carry 50% of the habitat capacity in different tints: a 
single toning (left map) and tinting of territories with different colors (right map) 

Note that the maps in Figure 7 and Figure 8 differ rather significantly. This can be 
explained as follows. The density of the Earth's population, which was used as the weight in 
formulas (14), varied in the range [1,04;10^]. The density of the capacity of the habitat, 
which was used as the weight in formula (15), varied in the range [0,77;15,26]. A couple of 
bands differ from each other in 690 times. In other words, the population of the Earth due 
to high density, i.e. by concentrating in large urban agglomerations, "economize" on space. 
Furthermore were calculated square accommodations 50% of the population, ^population 
and square seats concentration 50% capacity, ^capacity' represented in Figure 7 and 8. It 
turned out that 5’popuiation = 0,5131/?^, ^capacity = 2,1120 R^, where R is the radius of the 
Earth's. Thus, the surface of living 50% of the population in 4 times less the surface of the 
positioning 50% of the habitat capacity. 


Maps in Figure 8 can be viewed from a different perspective by constructing a density 
gradient field of the habitat capacity. We define in the spherical coordinate system a two- 
component gradient vector g = Vp of the density of the habitat capacity in the form: 


g = 

ocp cos (p OA 


(16) 
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The density of the habitat capacity is considered on a grid with a resolution of half a degree. 
To approximate the derivatives in (16), we take the four-point pattern, then we write 

(^). , “ ^(Pi+ij+Pi+2j—Pi-ij—Pi-2j)» (17) 

where i = 3, ...,357; j = 1, ...,720 and 

where i = 1, ...,360; j = 1, ...,720, it is assumed that in (18) a periodic extension with 
respect to the index/with period 720. 


Figure 9 shows the gradient field as a set of arrows, gij = of fh® density of 

the medium, calculated by formulas (16) — (18) with density of the habitat capacity (3) 
with parameters from point #5 of Table 1. 

In Figure 9, clearly in the form of a "fringe” of arrows, a noticeable increase the density of 
the capacity habitat is seen as it approaches to the shoreline, which once again confirms the 
well-known fact that half of the world's population lives in the vicinity of the coast of the 
world ocean. 

The field density gradient of the capacity of the habitat 
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Figure 9. The field density gradient of the capacity of the habitat 
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SECTION IV. 

CUTTING TERRITORIES IN TERMS OF “(not)HIGH - (not)FAVORABLY" 


The studies presented in the previous section allow the classification, cutting of territories 
in terms of "high - low", "favorable - unfavorable”. The term "high-low" means literally the 
following: the locality of residence is at an altitude higher than a certain threshold value Zp 
or, accordingly, less. The term "favorably - unfavorable” means that the density of the 
habitat capacity exceeds a certain threshold value Pp or, respectively, less. 

According to the previous section, the threshold values Zp and Pp were found from the fact 
that the capacity of the habitat of the territories covered by these conditions was 50%. 
According to Table 5, Zp = 820 m and Pp = 1,82. 


Consider now the conditional space with the coordinates p (0 < p < p^ax) Z (0 < Z < 
Zmax) draw two lines p = Pp and Z = Zp, then any area of habitat on the Earth can be 
assigned to one of the four areas on plane with coordinates pand Z. Under coordinate p we 
understand the density capacity of the habitat, and a coordinate Z — height of the location 
above sea level. Figure 10 shows the positioning of each of the four areas with an 
approximate preservation of the scales by coordinates. 



IT I 50%!/ 




Pp = 1,82 


^max 


= 8,39 


=> 

P 


Figure 10. Positioning of each of the four areas in the space with coordinates p, Z 


According to Figure 10, the area I is determined by a pair of conditions: p > Pp and Z < Zp; 
area II — p <Pp and Z < Zp; area III — p < Pp and Z > Zp; area IV — p > Pp and Z > Zp. 
In terms of "highly favorable”: the first area is most conducive to living, i.e. these territories 
are not very high, and they are climate-friendly; the second area corresponds to areas that 
are not very high, but they are unfavorable from a climatic point of view; the third area 
corresponds to the most severe territories, located high and which are unfavorable from 
the climate point of view, finally, the territories of the fourth territories, although they are 
located high, are favorable from the climate point of view. Figure 10 also shows the 
coloration specific for each type of territories used in constructing the map in Fig.l2. 
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Figure 11 shows four maps of the territories corresponding to each of the four squares in 
Figure 10. Territories within a separate category are colored in Figure 11 by random 
colors. The initial data were taken according to calculations of the previous section. The 
results of calculations of the areas of territories, volumes of capacities of the environment 
and the number of territories for each of the types are given in Table 6. 



Area II (low and unfavorable) 

100 r 



100 I-,-,-,-.-.-,-,-, 

-200 -ISO -100 -50 0 so 100 ISO 200 



Area I (low and favorable) 



Figure 11. Maps of the territory with parameters of density of capacity of habitat and 

height of positioning according to four areas 

The squares of areas in the Table 6 are given in units of R^, where R is the radius of the 
Earth. In particular, if the whole territory of the first territory is selected as the unit, the 
ratio of the areas of all four types of territories will be 1: 0,59: 0,35: 0,30. Finally, if the total 
capacity of the habitat of the territories of the first area is chosen as the unit, the ratio of 
the capacities of the habitat of all four types will be 1: 0,41: 0,21: 0,38. 


Table 6. A set of areas (5), volumes of habitats (f/) and the number of 
_ territories (IV) for each of the areas _ 


Variable 

Areal 

Area 11 

Area III 

Area IV 

Sum 

S/R^ 

1,7909 

1,0632 

0,6276 

0,5420 

4,0237 

U 

3,8525 

1,5796 

0,8095 

1,4623 

7,7039 

N 

391 

473 

401 

454 

1719 


In Figure 12,a shows a single map for the territories of all four qualities, while the site of 
the first area are colored green, the second area is blue, the third area in violet and finally 
the fourth is yellow. 

Let us consider in more detail the map depicting the state border of the Russian Federation 
(for 2008), as well as the territories classified according to Figure 10. The result is shown 
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in Figure 12,h. It is clearly seen from Figure 12,h that practically the entire territory of the 
Russian Federation is painted in blue, it is classified as "low - unfavorable”. A small 
fragment of the territories of the first category "low - favorably” (painted in green) is 
located in the west of the country. 



Map of territories of all four qualities Hith marking of the state border 
of the Russian Federation (for 2(M)8) 


Map of the positioning of the territories of all four areas 



Figure 12,a. Map of joint positioning of 
territories of all four areas 


Figure 12,b. Map of territories of all four 
qualities with marking of the state border 
of the Russian Federation (for 2008) 


We will calculate the so-called "index of the diversity of territories”, D for all countries and 
territories, the list of which was considered earlier in section III. To do this, we define the 
concept of "territory index”, h, which can take values from the set (1; 0,41; 0,21; 0,38), 
introduced above in connection with the mutual subordination of the habitats of four types 
of territories. Define the index of the territory for each of the points of the grid 360x720, 
then we obtain the matrix hij, i = 1, ...,360; j = 1, ...,720. In Section III, the database for 
2008 was reviewed for 246 states and territories for which we will find the indices of the 
diversity of territories D^, k = 1, ...,246. 

The index of diversity in an arbitrary territory is defined as the variance of the values of the 
territory index, i.e. 


Dfc = var{hi j, ((pi,Aj) G i = 1 ,... 360,j = 1, ...,720), (19) 

where var — variance operation, — area occupied by the k-th state or territory on the 
Earth sphere, k = 1, ...,246. 

The territory diversity index (19) was calculated for all 246 States and territories. Figure 
13 shows a fragment of the summary data for the selected set of States. If, for example, to 
compare the indices of the diversity of the territories of the USA and of the Russian 
Federation, this ratio will be equal to Dusa/^rf = 4,894. Thus, the territory of the Russian 
Federation is about five times less diversity then the territory of United States. This means 
that the management of the US territory as an integrity is more complex than the similar 
management of the territory of the Russian Federation. 

As a result of the studies presented in this and the previous sections, the following was 
revealed. Most of the population of the Earth seeks to live in territories classified as "not 
highly - favorable” or, in other words, closer to the land-ocean interface. This half of the 
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population is the most socially mobile, because it is oriented to the ocean, more precisely, 
to ocean traffic. Another, more conservative part of the population, is oriented to land and, 
consequently, continental traffic. These considerations reproduce the opposition of two 
political doctrines conventionally called "thalassocracy” and "tellurocracy”, well known in 
geopolitics (H.J. Mackinder [7], K. Haushofer [8] and several other authors), the first of 
which is oriented toward acquiring political power through the ocean, and the second — 
through the continent. 


The values of the diversity index of the territories of the selected set of States 

■k United States 
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k Iran (Islamic Republic of) 
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k Japan 
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★ Germany 
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Numbers of States and territories 


Figure 13. The values of the diversity index of the territories of the selected set of 

States 


In this section, the cutting of land territories in terms of "(not)high — (not)favorable" is 
presented (Figure 12). This cutting can be considered as a prerequisite or a framework for 
explaining the existing configuration of State boundaries. It is clear that in the final 
explanation of the configuration of States borders there is always a fair share of the 
prehistory of the formation of a State (territory) in a single context of the world political 
system. In this section, we were also interested in the unification and classification of 
territories in terms of climatic data and relief 


Consider in the future a different method of cutting of territories, in which the decisive role 
is played by traffic. This way of regionalization should take into account both geopolitical 
features (thalassocracy and tellurocracy) of territories management. 

In the beginning, we will develop an algorithm for the random distribution of a set of points 
on the land surface, taking into account the density of the habitat capacity. In the future, we 
will determine the traffic between points, the cost of traffic and, accordingly, determine the 
task of minimizing the cost of traffic by optimally moving the points serving as logistic 
nodes between the territories. The optimal position of the set of communication nodes and 
the coastline will determine the desired territorial subunits — geopatoms. 
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SECTION V. 

RANDOM DISTRIBUTION OF POINTS ON THE EARTH'S SURFACE 


We represent the integral in (8) as a sum in the following way. Define within the 
environment on the surface of the Earth N points with the coordinates Zj = (q)i, Ai), i = 
We will assume that all points are localized within the land. Let each point Zj 
corresponds to some region (Ot- As a result of the above considerations, the global capacity 
of the habitat (8) can be rewritten as: 

U = 'Zi^^Ui,Ui = ff^ dcpdAcos(p ■ p. (20) 

Using the Monte Carlo method, we will play the procedure of random application of N 
points within the land. It is necessary to play a two-dimensional random variable z = 
(^,/I), containing the components of latitude and longitude. Taking into account the 
previous sections, consider a grid of latitude and longitude values with a resolution of 0,5®. 
In this case, the two-dimensional random variable z will take the set of values z^j = 
i = 1, ...,360; j = 1, ...,720. We write the distribution matrix Pij of the random 
variable z in the form: 


Pij =p{z = 



Pl,j cos (pi 


V360 

^k=l 


ti=lPk,lCOS(pk 


( 21 ) 


where i = 1, ...,360; j = 1, ...,720. 


If we sum the matrix (21) by rows and columns, we find the series of latitude distribution, 
Pfp i and longitude, pxj, i.e. 


V<p.i = Y/ji\ Pij ,i = l .360; pxj = 1.^1 Pij .j = 1 .720. (22) 


We construct series of distributions (22). To do this, we substitute in (21) the density of 
the capacity of the habitat according to formula (3) with the parameters from point #5 of 
Table 1. The result is shown in Figure 14. 

The study of the left-hand graph in Figure 14 confirms the expectation that the probability 
to the poles should decrease with the random play of the latitude component of the random 
vector z. Taking into account the form of the density of the medium (3), it is clear that the 
components of the random vector z are not independent. For this reason, in order to play 
each of the components, it is necessary to know one of the two conditional probability 
matrices of the form: 


_ Vi,j _ Pij 

P<P\Uj - 


(23) 


where i = 1, ...,360; j = 1, ...,720. 
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Figure 14. The series of distribution of latitude components (left graph) and 
longitude (right graph) of the random two-dimensional vector z 


Consider the procedure for determining each of the components of the vector z = {cp, A) 
individually. Taking into account (22), the random variables ^and A assume the values cpi 
and Aj with probabilities p^pi and pxj, where i = 1, ...,360; j = 1, ...,720. Prepare sets of 
cumulative probabilities of the following forms: 

Pzcp.i = Sfc=i V(p,k , t = 1, - ,360; p^xj = PA,k J = 1, - ,720. (24) 


Taking into account (24), determination of the value of any of two random discrete 
quantities (p, A reduces to the following procedure. For definiteness, we choose the latitude, 
for longitude the procedure is similar. We generate a random number ^ uniformly 
distributed on the interval [0;1]; find the number of the interval from the set [0;px<p,i], 
[Pi.(p,i. Pi.<p. 2 \, -, [Pz<p. 359 > Pi:(p. 36 ol in which a random number falls. Denote the resulting 
number by the expression i(^), then the played value of the random variable ^ will be (Pi(^y 


Let's go back to determination of the values a two-dimensional random variable z = (<p, A). 
Suppose, for example, that the latitude is determined out and that the value is 
obtained. We need to determine the longitude at a certain latitude i.e. it is necessary 
to use the conditional probability Px](p,i.j from (23), when i = i(Q. In other words, the 
problem reduces to determining a conditional discrete random variable A\(p = that 

takes the values Ai, ...,^720 with probabilities Px\<p,i,iiO’ PM<p.72o.m- ^ 

procedure similar to the latitude determination, i.e. find the corresponding cumulative sets, 
generate a random number rj uniformly distributed on the interval [ 0 ; 1 ] and determine the 
interval number y (p), which gets a random number rj. 

As a result of the above algorithm, the procedure for determining the two-dimensional 
random variable z = (q),A) is completed, and generation of the pair (<Pj(f), Ay(^)) can be 
continued without restrictions. In specific calculations, the algorithm was somewhat 
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different. Let N value pairs ^ ~ be generated. For even values of n, 

the latitude was first determined, and then, using the conditional probability, longitude was 
determined and, conversely, for odd n, longitude was determined first, and then with the 
conditional probability the latitude was determined. 

Figure 15 shows the result of the application described above in the form of one of the 
versions of the random application of iV = 10^ points within the land, taking into account 
the density of the habitat. 

The random distribution of points, taking into account 



Figure 15. The random distribution of iV = 10^ points, taking into account the 

density of the capacity of the habitat 
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SECTION VI. 

TRAFFIC BETWEEN TERRITORIES 


Consider the presence of correspondence or traffic between a pair of territories with 
numbers i and j {i,j = 1,..., N). In the previous section, an algorithm was developed for the 
random distribution of N points on the Earth's surface, taking into account the density of 
the capacity of the habitat. With each of the points we associate a certain territory, which 
will have corresponding a capacity habitat Uj, i = As such territories, we take 

Voronoi polygons for a finite set of points i = on the surface of the 

Earth. A set of N points can be interpreted, among other things, as logistics nodes for 
transporting part of the habitat capacity between territories. 

Figure 16 shows an example of the distribution by the Monte Carlo method of Ai = 10^ 
points on the Earth's surface, taking into account the procedure (21) — (24). There, along 
with the shoreline and poles, the corresponding Voronoi polygons are plotted. When 
constructing the Voronoi polygons in the latitude-longitude coordinates, a singularity 
associated with the poles took place. Four additional "fictitious" points were introduced: 
Zn+1 = J(900,180°), z ^+2 = /(-900,180°), z ^+3 = /(-900, -180°), z ^+4 = 

/(90°, —180*^), where x = const > 0. In the calculations, the parameter was chosen equal 
to 5. Next, Voronoi polygons were found from the extended set of points = (<^j,Aj), i = 
1, ...,N,N + 1,N + 2,N + 3,N + 4. In the framework of this procedure, Voronoi finite 
polygons are guaranteed to include the coastline, and unbounded polygons were moved to 
the fictitious area. 


The map of the Earth’s surface with a random collection of points within the land and Voronoi polygons 




-60 
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Figure 16. The map of the Earth's surface with a random collection of IV = 10^ points 

within the land and Voronoi polygons 
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In the beginning we will define the traffic, escort of which is cost-free. We define the matrix 
ai j, i,j = 1,..., N, the elements of which characterize the transfer of the i-th fraction of the 
habitat capacity to the ;-th fraction of the habitat capacity in connection with traffic. The 
matrix i,j = 1, ...,N acts as a measure for the exchange of the capacity of the habitat 
between a set of points Zj = ((Pi,A.i), i = 1, ...,N on the Earth's surface. According to the 
definition, it follows that 

= 1, (25) 

where an — the share of the i-th capacity of the habitat is not subject to traffic, i.e. this is 
the fraction that remains inside the i-th territory. The set of capacities of the habitat Uj, i = 
1, ...,iV acts as some resource that in the context of traffic is considered universal and 
additive. 


We find the amount of habitat capacity, [/•*' ^ left due to traffic the i-th territory and the 
corresponding amount of habitat capacity, arrived in the i-th territory, then 






jjW _ yN 

' ^j= 


lj:^i 


(26) 


Taking into account (25), it can be established by direct verification that there is a balance 

of arrival and departure of the habitat capacity, 

i.e. due to the cost-free traffic, how much it takes so much and comes. 


Now, if we move a part of the habitat from point i to point j, the share of the resource yij 
due to traffic is lost, we assume that 0 < Yij < 1, i j] i,j = 1, —, N. We believe that within 
the logistic point there are no costs, i.e. Yu = 0, i = 1,..., N. The matrix Yij, i,j = 1. —, N is 
called the transport cost matrix. Taking into account (26) we will form the functional, Tr of 
all transport costs: 






N 




(27) 


The meaning of the transport functional (27) is transparent: when the resource is moved 
from point i to point j, its part in quantity Yij is consumed, similarly, when the resource is 
moved from point j to point i, the resource part in quantity yj i is also expended. The 
transport functional (27) can be simplified somewhat, taking into account that the first and 
the second double sums are the same, this is seen after renaming the indices, i.e. 

Tr = 2 oCijYijUi- (27') 

Using the formula (27') we can put the optimization problem of minimizing the cost of the 
traffic 7’r(zi, ...,z^) by a suitable displacement of the z^, ...,z^ points within the land. It is 
clear that the global minimum cost of transport is zero value, when Tr = 0. In a 
mathematical form, the minimization problem reduces to finding a set of points {z 1 , 
which may not be unique, in which one of the local cost minimums for traffic is realized, i.e. 
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( 28 ) 


Tr(zl = min rr(zi.z„). 

zi,...,zjv 

For further progress in solving of the optimization problem (28), additional clarification of 
the form of the matrices ai j, Yij, i,j = 1, , IV is needed. Let us consider the well-known in 
the theory of transport flows the so-called gravitational model [9, 10], in which some 
generalized "distance" di j (dij > 0) is introduced between the i-th and ;-th territories. 
Note that the distance matrix djj, i,j = 1, in general, is not symmetric, i.e. dij ^ dj i, 
whereas by definition it is assumed that dj ; = 0, i = 1, Note that the condition d^j = 
0 ,i ^ j does not necessarily imply that the and Zj points are the same in the physical 
space. 

We choose the exponential dependence of the coefficients of the matrix i,j = 1, 
on the distance, then we can write the following representation: 

- yN g-pdiy ( 2 ^) 

where yd — some nonnegative parameter. One can verify by direct verification that the 
condition (25) holds for the matrix i,j = 1,..., N in the form (29). 

Following the exponential dependence of the matrix a^j, i,j = 1, ...,N on the distance in 
(29), we choose the following simplest representation for the transport cost matrix: 

Yu = ir(l - (30) 

where r and P 2 — some constant nonnegative coefficients. According to (30), for djj = 0 
there is no transportation cost, i.e. Yij = 0. Finally, for ^2 ^ and for dij 9 ^: 0 it follows 
that Yij ^ i.e. the share of transport costs becomes a constant value of Ir. The last 
remark implies that 0 < r < 2 . 

We substitute (29), (30) into (27'), then we find the following expression for the total 
transport costs: 

Tr = ri:f=i[l - (27") 

The transport costs in the form (27") are characterized by the following properties, which 
are checked directly. First, when the distance between the points becomes zero, i.e. d^j = 
O.i.j = there are no transportation costs, Tr = 0. Secondly, when the distance 

between points tends to infinity, i.e. dij 00 , i ^ j; i,j = 1,..., N, the transport costs due to 
the selected dependencies (29), (30) also tend to zero, Tr 0. Note that a distance of zero 
or infinity between two points does not mean that the points merge or diverge to infinity. 
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Consider transport costs in the form (27") as a function of the parameter ^ 2 < i-e- Tr = 
Tr^Pi). In this case it is obvious that rr(0) = 0. Now let ^2 ^ then rr(^ 2 ) ^ Tr(co) = 

For the last transport cost functional, it is obvious that the 
global minimum equal to zero is achieved only when the distances between points tend to 
infinity. The latter case will be considered in the future from the point of view of 
computational experiment. In this case, the points cannot come together, they will repulse 
and fill the maximum of the habitat, because they are "locked" on the surface of the 
terrestrial sphere. 

The condition of the transition limit ^§2 ^ we call "the minimax transport doctrine”, 
which means the minimum transport costs at a maximum filling capacity of the habitat. 

We substitute (27”) in (28), then in the general case the problem of minimizing transport 
costs will be rewritten in the form: 


Tr(zl....,z*^)= min 


N 


[1 


yN , 




yN 

^k=i^ 


-P^i,k 


]u,. 


(31) 


To clarify the problem of minimizing (31) in a simplified form, when it is assumed that 
1^2 00, i.e. there is a minimax transport doctrine, then 


Tr(zl,...,z*^)= min r2f=i(l- 








)Wr (31') 


In order to determine the algorithm for calculating the distance matrix di j, i,j = 1,..., N we 
will give a series of physical considerations on the energy costs of moving one unit of cargo 
weight from = (?>i, A^) to Z 2 = {(P 2 ’^ 2 ^- Let the route of movement of freight is defined 
in the form of a line: cp = <p(s), A = A(s), 0 < s < 1. The argument of the line 
parametrization is the length of the line, s, changing from zero to its maximum value, I, 
equal to the length of the route, while it is assumed that cpi = (p(0), A^ = A(0) and ^2 = 
(p(l),A2=A(l). 


Let in the beginning the route of movement is completely located on the surface of the land. 
In this case, the energy cost of moving one conventional unit of cargo weight consists of 
three characteristic variants of movements: 1) vertical movement of cargo upwards due to 
the peculiarities of the relief; 2) vertical lowering of cargo due to the peculiarities of the 
relief; 3) movement of cargo along an inclined surface. The first two points of the variants 
of the movement describe the energy costs of moving the load "up and down". The last 
point is mainly characterized by energy costs to overcome rolling friction within the 
framework of such modes of transport as road and rail. Let Z = Z(^, A) be the relief of the 
land surface, then for the selected route of movement of the conventional unit of cargo 
mass, we can write the function Z(s) = Z(^(s), A(s)), 0 < s < 1. 


Taking into account the stated physical considerations, we will write the formula for 
calculating the distance along the chosen trajectory between a pair of points z^ and Z 2 , 

dearth- 
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dearth = 91 /q 92 /q + ^3 /q “ (f)^ds, (32) 


where K(t) = 1, t > 0; K:(t) = 0, t < 0 — so-called "single" function. The unknown non¬ 
negative parameters 9 i, 92>93 characterize the contribution of each of the types of 
movements for the displacement of one conventional unit of cargo mass. The integrals in 
(32) are called transport integrals. 

Let now the route of movement of one conventional unit of mass of cargo lies entirely in the 
sea, i.e. the points of departure, arrival, and all other points of the route lie on the surface of 
the sea. In this case, the energy costs associated with the movement of one unit weight of 
the conventional cargo by sea transport are associated with the overcoming of viscous 
friction. To calculate the distance along the selected path between a pair of points and Z 2 , 
dsea oiie can use formula (32). We believe that water transport moves along a horizontal 
surface, for which it can be assumed that dZ/ds = 0. Taking into account the last integral 
in (32) and performing elementary integration, we write the corresponding transport 
integral in the form: 


^sea “ 9 ^ 1 , (33) 

where — a certain nonnegative parameter that takes into account the averaged features 
of viscous friction in the water environment of the aggregate water transport. 

Note that an arbitrary route between points of departure and of destination can be broken 
down into stages of the movement only on land or only at sea. Applying either formula (32) 
or formula (33) to each step and adding the values obtained, we find the total distance 
between a pair of points. Distances calculated by formulas (32), (33), are not distances in 
the usual sense of the word. They rather act as effective distances, which can always be 
measured by calculating the average energy of the movement of one conventional unit of 
mass of cargo from the point of departure to the destination. 

Based on physical considerations, we estimate the order of magnitude of the parameters 
9i>92>93>94- When lifting a load on an inclined plane, work in the part of lifting is 
expressed by the formula: mg^ds, where m — cargo mass, and g — acceleration of 

gravity. Let us find the specific work per unit of cargo weight by dividing by mg, it is ^ ds, 
i.e. g^ = 1. 

A somewhat different situation when moving cargo along the inclined plane in the case of 
descent. It is not clear how much the average energy for conventional land transport (road 
and rail) is required for descent. The simplest case is to assume that the same amount of 
energy is required as in the ascent, then we find g 2 = 1- 

To evaluate the third parameter, it is necessary to take into account the formula for 
determining the friction of rolling. Let / be the friction of rolling, and R is the radius of the 
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wheel, then we can write the following formula for working on overcoming the rolling 
force: ^mgji-(^yds. Note that the friction of rolling rubber on asphalt for motor vehicles, 

/auto the friction of rolling steel on steel for rail transport, /ra,] can be considered equal 
in order of magnitude, i.e. /auto — /rail — mm. Consider that the average wheel radius 
total car and wagon railway transport order 0,5 m. In this case, the specific work on moving 

the unit of cargo weight will be Substituting / = 0,5 ■ 10 and R = 0,5m, we 

find ^3 = 10 “^. 

We now proceed to estimate the numerical coefficient ^ 4 . To do this, calculate the force of 
resistance R for a typical container ship [11] with dimensions: length, L = 157m; width, 
B = 25,4m; draught, T = 9,22m, displacement, D = 22800 ton; design speed, v = 

19 knots = 9,77m/s. To calculate the resistance force of friction it is common to use the 

2 

formula R = ^^11, where ^ is the complex aggregate friction coefficient, p is the seawater 
density, fi is a wetted surface of ship's hull. We divide the force of resistance by the weight 
of the vessel, equal to mg = p^Vg, is an average mass density of a ship with a cargo, V is 
the volume of the ship, approximately its displacement, g means the acceleration of gravity. 
As a result, we can write the following estimate for the sought coefficient: g^ = We 

believe that ^ = 10“^, p~p*. The value of the wetted surface is calculated by the formula of 
Muragin D. = LT ^1,36 + 1,13 x 0,62 x — 4,76 ■ lO^m^. As a result, we find g^ = 10“^. 

We note that within the framework of our analysis, in order of magnitude, the specific 
energy costs for the movement of a unit of cargo mass with the help of land and water 
transports are the same. 


In the optimization problem (31') a set of capacities of the environment of individual 
regions Uj, i = were considered as fixed. However, this does not correspond to the 

procedure of cutting areas using Voronoi polygons. Indeed, with the optimization of the 
location of the points z^, ...,z^, their positions will change, which will lead to a change in 
the arrangement of the corresponding Voronoi polygons. For this reason, in (31') it is 
necessary to keep in mind the set of Uj, i = 1,..., N, obtained in connection with the current 
set of points z^, ...,Zyv and the fact that Y.i=iUi = U, where U is the capacity of the habitat 
within the land, Ui is the capacity of the habitat of the Voronoi polygon a)^. As a result, the 
optimization problem (31') should be rewritten as: 


Tr(zl, ..., 2 *^]) = min r2f=i(l- 




U = Ui,Ui = ff dcpdX cos (p ■ p. 


)«i. (34) 

(35) 


To solve the optimization problem (34), (35), we will use the simple iteration method, 
choosing the set Ui, i = in (34) from the previous iteration. The search for a 

minimum in (34) allows us to find the set of points z^, ...,z^ at the next iteration, by which 
we define a new set Ui,i = 1, ...,N and so on. Let rbe the iteration number, t = 1,2,..., then 
problem (34), (35) can be rewritten as: 
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Tr(z[^\ = minrSf=i(l- 


,(t-1) ,,(t-1) _ 




k=i 




(T-1) 


(34') 


U = ,w,- — ff (T-i)d(pdAcos(p ■ p. (35') 


The optimization problem (34'), (35') is solved iteratively until the sequences 
z-^\z^^\ ...,z[^\... converge for each i = 1,...,N. It is assumed that the set 
determines the initial location of the points, and the collection uf^\ i = 1,...,N is the 
corresponding list of capacities of the Voronoi polygon environment. 
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SECTION VII. 

THE METHOD OF CALCULATING OF THE DISTANCE MATRIX 


To calculate the generalized distance (32), (33) between an arbitrary pair of points on the 
surface of the Earth = ((p^, X{), = ijp^, A 2 ) construct in the beginning the trajectory 

of a small fragment of the total circle of the Earth passing through a given pair of points. We 
assume that the total circle lies in the plane passing through the center of the Earth. In 
other words, we construct a line: (p = (p{s),X = A(s), where s is a parameter describing the 
length of the fragment of the circle. We assume that 0 < s < 0, ^(0) = ^i,A(0) = and 
(p{d') = (p 2 ,X{d') = X 2 , where 6 is the arc length of the fragment of the total circle of the 
Earth. 

We define two unit vectors of length n^, n 2 , which point to a pair of selected points = 
((Pi,Ai),Z 2 = (^ 2 . ^ 2 )-then 


rifc = (cos Afc cos (pk , sin cos cp ^, sin cp^'), k = 1,2. (36) 

Let the vector n = (Ux, Uy, n^) of unit length indicate an arbitrary point of the complete 
circle passing through the selected pair of points. It is clear that the vector n lies in the 
plane formed by the vectors ni,n 2 . Taking into account (36), after simple transformations, 
we find 


n = n(s) 


sin(0-s) 
sin G 


ni + ^n2, 

sin0 


(37) 


where ^is the angle between a pair of vectors n^,n 2 . According to (37), it is obvious that 
n(0) = rii and n(0) = n 2 . 


Taking into account (37), and also considering that n(s) = (cos A(s) cos (^(s), 
sin A(s) cos (p{s ), sin (pis)), we find the parametric recording of the complete circle in the 
coordinates "latitude - longitude”: 


(p = (pis) = arcsin(nz) = arcsin[ ^‘"[^g^^ sin + 


sin 6 


A = A(s) = arctg ^ = arctg 


sin ^ 2 ]- 

sin( 0 - 5 ) sin Ai cos (pi+sins sin/l2 cos (p2 
^ sin( 6 -s) cos cos <^i+sins cos A2 cos <^2' 


(38) 

(38') 


Note that the formula (38') is true when % > 0. In the other two cases: 1) % < O.riy > 0; 
2) rix < 0,ny < 0 to the expression in (38') it is necessary to add and subtract n, 
respectively. 

When carrying out a line between a pair of points, according to formulas (38), (38'), it is 
also necessary to find its intersection with the coast line to identify those parts of the 
trajectory that lie separately on land and at sea. To depict the trajectories (38), (38'), two 
cases must be distinguished: 1) |Ai — A 2 I < tt; 2 ) |Ai — A 2 I > n. In the first case, the zero 
meridian will be placed in the center of the map, i.e. the Atlantic Ocean will be located in the 
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center. In the second case, after shifting the longitude A A + n, we construct a map with 
the center along the line of the date change, i.e., with the center in the Pacific Ocean. 

Figure 17 shows an example of the positioning of ten points on the land surface. The points 
are chosen randomly according to the algorithm of section V, i.e. taking into account the 
density of the capacity of the habitat. The points are numbered and are shown in each of 
the two cards. In Figure 17, taking into account the spherical geometry, all 10x9/2 = 45 
lines are constructed, connecting each pair of points. In addition, the lines are marked in 
part of their passage overland (solid line), and by sea (dashed line). According to Figure 17, 
some of the binary lines are positioned on the right map with the center in the Pacific 
Ocean. 


A set of small fragments of complete circles A set of small fragments of complete circles 




Figure 17. A set of small fragments of complete circles connecting pairs of points 

According to formula (32), to calculate the effective distance on the displacement of one 
conventional unit of the cargo mass on the Earth's surface, it is important to know the 
derivative of the relief Z = Z(q),A) along the route, which we describe with some 
trajectory. Let, for example, a path be established between a pair of points = (^i,Ai), 
^2 = such that (p = (p{s),A = A(s), where s is a parameter describing the distance 

traveled from the starting point. In this case it is obvious that 

As a route between a pair of points on the surface of the Earth's sphere, we choose a 
smaller fragment of the total circle passing through a pair of points. A suitable trajectory is 
presented in the form of formulas (38), (38'). It remains to find the derivatives. After 
simple calculations, we get: 


d(p _ - cos(0-5) sin</?i+cos5sin</?2 

ds ^sin^ 6-[sin(6-s) sin<^i+sin5 sin 
dX _ sin 6 sin(/l 2 -Ai) cos cp^ cos (p 2 

ds sin^ 6-[sin(6-s) sin</?i+sins sin (pzV 


(39) 

(39') 
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Let us describe a large fragment of a complete circle passing through some pair of points on 
the Earth's surface. Taking into account the description of the small fragment of the 
complete circle (37), (38), (38'), it is expedient to parameterize the large fragment in two 
steps: 1) for values of the arc length parameter s from the range [jU, 0], where jU < 0; 2) for 
the values of the arc length parameter s from the range [6,211 + fi]. We find the angle fi 
from the condition that the longitude (38') coincides with the meridian of the date change, 
i.e. A(ju) = —n. From the solution of the last equation, we find n = iiq, when (Uq < 0 and n = 
Hq — n, when jUq > 0 , in this case 


(Uq = arcsm 


sin d sinAi coscpi 




Figure 18 shows an example of a complete circle passing through a pair of randomly 
selected points on the Earth's surface. The circle is divided into three arcs: 1) a small 
fragment of the full circle is marked with markers in the form of stars; 2) a large fragment 
of the total circle parametrized by a segment of the arc s from the range [jU, 0], where n < 0, 
is marked with pentagrams; 3) a large fragment of the complete circle, parametrized by a 
segment of the arc s from the range [6,2n + ju] and is marked with hexagrams. Each of the 
arcs is marked with an appropriate angle. The arrow that separates the pentagrams and 
hexagrams indicates the date change line. After combining the second and third fragments 
of the arcs, we obtain a large fragment of the complete circle. 


Partitioning the entire circle into three fragments 



Figure 18. Partitioning the entire circle into three fragments 

We introduce the distance matrix D = {di j,i,j = between all pairs of points. 

Define similar distance matrices found for small Dgmau = {dsmaiuj’^J = —,N} and large 
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D\arge — {d]arge,ij’^J — fragments of arcs of corresponding complete circles. 

Considering the solution of the problem of minimizing the cost of traffic, assume that the 
required distance matrix D is the element-wise minimum of the pair of distance matrices 
^smaii and i.e. 

D lTlin(f^small'^large) ~ rnfr(t^small,ij'^large.ij) » ~ 

We perform a computational experiment to estimate the mean value of transport integrals 
in formulas (32), (33). To calculate the transport integrals, we need the partial derivatives 

d Z d Z 

— and —, which we compute using the finite differences in the form (17), (18) with the 
resolution of the relief 0,5°; as well as ordinary derivatives — and —, which we find 

ds ds 

according to (39), (39'). In the rest the fragment of the full circle, connecting a pair of 
points, was divided into parts, passing separately by land and by sea. For each part, a finite 
difference grid was constructed with the number of points equal to the value [2 -I- 

0,25^ (^2 “ (PiY + (.^2 ~ where [... ] is the function of the integer part of the number, 
and the pairs (^i,Ai) and ((P 2 ,^ 2 ) denote the beginning and the end (in degrees) 
considered parts. In the future, in some cases, in order to save computing resources, 
parameter 0,25 was replaced by 0,125. As a result, for each pair of points from the total 
number N x (N — l)/2, we compute transport integrals and find the mean values from 
each set. 

Table 7 summarizes the computational experiment in a format where, randomly, according 
to the algorithm of Section V, N = 320 points were selected and the matrix £>smaii was 
calculated. The top line in the first row of Table 7 denotes the averaging operation. Taking 
into account the estimates from the preceding section, fr follows 

that each term in (32), (33) makes a comparable contribution from the point of view of the 
mean values. As a result, let us consider the choice of the values of the parameters = 
92 = '^.93 = 10 “^, .04 = 10 “^- 


Table 7. Average values of transport integrals 


1 

1 

1 


I dZ AZ., 

J 

\ dZ ( dZ\, 

-J 


1 

0 

0 

0 


7,26 ■ 10-4 

4,85 ■ 10-4 

0,7556 

0,6383 


According to the definition in (32), (33), and in Table 7, the integrals are taken along a path 
connecting a pair of points on the surface of the Earth. Since all the routes are located on 
one of the fragments of the complete circles, the transport integrals can be approximately 
considered proportional to the length of the arc. From this point of view, on the average, we 
can assume that the distance matrix D is the same as the distance matrix found for small 
fragments of complete circles, i.e. we assume further that D = Z)smaii- 
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SECTION VIII. 

THE TRANSPORT COST MINIMIZATION ALGORITHM 


Let us return to the problem of transport costs optimizing in the form (34'), taking into 
account: 1) the features of the "land - sea” geometry and 2) the computational complexity 
of the transport cost functional Tr calculating and, in particular, its partial derivatives of 
variables z^, Suppose that in the previous step of the iterative procedure, according 
to formulas (35') the volumes of the habitat capacities are calculated, which allows to 
determine the transport costs functional (34') within the current iteration. 

To find the minimum of the function Tr = Tr{(p^, Ayv) we consider a slightly 

modified gradient descent scheme. Let us write the following system of ordinary 
differential equations: 


dcpi 

dt 


-hySign(S,^ = -aysign(^),i = 1. N, 


dcpi dt 




(40) 


where t — some auxiliary argument, sign(x) = 1, x > 0; sign(O) = 0; sign(x) = — 1, x < 0 
— so-called "sign” function, and the sets of non-negative coefficients {%, ...,ayv} and 
(hi,..., hyv) are defined below. 


We replace the ordinary derivatives in (40) by finite differences, then 


6 (pi = -5t ■ hiSign(^), dAy = -5t ■ aysign(^), (41) 


where Sq)i,SXi,i = 1,... ,N are the finite increments of functions, and St is the final 
increment of the argument. 


Taking into account (41), we write the transition algorithm from the current values of the 
functions {(p^, ...,(pf^,X^, ...,Ayv} with the value of the argument t to the new values of the 
functions [(p [,..., (p'fj, X [,..., A^} for the value of the argument t -I- St, namely 

(p'i = (pi-St- hySign(^),A; = Ay - dt ■ aySign(^), (41') 


where i = 1, ...,1V. If after recalculation, according to the algorithm (41') one or more 
points are outside the land, we reduce the corresponding coefficients from the sets 
(ui, ...,ayv} and (hi, ...,hyv} so that the points return to the land. In case the coastline does 
not interfere with the movement of points, the corresponding coefficient values are 
assumed to be equal to one. 


The value of the parameter St in (41), (41') was chosen from the set of steps St^, —, of 
the form: 


Strn = Elq^,Stjn+:^ = ....St^ = "/q”, (42) 
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where S > 0, q > 1 — some constants, and m,n(m = 0, m<n) — natural numbers. 
Specific values of the parameters E, q and n were chosen in the course of numerical 
calculations. 

With allowance for (41'), (42), we calculate the new positions of the points 
X'k i, —, by the formulas: 

dTv dTv 

(pk,i = <Pi- ■ biSign(— ),Afc j =Xi- Stk ■ aiSign(— ), (43) 

with steps Stk, k = m, ..., n, as well as the functional values of transport costs: 

TVk '^^(xPk,i> ■■■' Vk,N’^k,i’ ■■■' ^ ■■■' (44) 

Assume that in the set (44), when viewed from left to right, there is a number k^, at which 
the transportation costs become smaller than those that took place in the previous step of 
the gradient descent scheme, i.e. 

TVk^ < rr(^i, ...,^yv.^i.= m,....n. (45) 

If inequality (45) is satisfied, for some /c* it is assumed that the gradient descent step is 
completed and the new positions of the points are considered equal: 

^1 — ^K.i’ —’V'n — —’^'n — ^kt,N- (46) 

If inequality (45) is not satisfied for all values of /c* = m ,..., n, then this step of the gradient 
descent procedure is also considered complete and the new positions of the points are 
chosen equal to: 


~ ^n,l’ ~ Vn,N’^l ~ ^n,l’ ~ ^n,N- ( 46 ') 


To close the minimization algorithm (40) — (46') it is necessary to clarify the procedure 
for counting the number m in the set (42). We denote the number /c* from the previous 
step of the procedure (41') — (46') with the symbol k^, then the number m in the set (42) 
is calculated according to the formula: m = k* = 0 and m = k* — 1, k* > 0. 


According to the algorithm (40) — (46'), it is necessary to have partial derivatives of 
transport cost functions Tr = rr((pi, on the coordinates of the points. 

Since it is difficult to do this analytically, it makes sense to consider the finite differences. 
We introduce some increments of the arguments A^^AAj,i = 1, ...,iV and write the 
following expressions: 


Pi = Tri... , (pi + A(pi, ..., X), Ri = Tr((p, ..., Aj + AAj,...), (47) 
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where i = 1, In this case, as above, we will assume that N points with coordinates 
(jPi + A(pi,Ai + AAj), i = are located within the land. If this is not the case, we vary 

the increments /\(pi, AAj, i = 1, Taking into account (47), we assume that 

dTr Pi-Tr dTr Rt-Tr 

7 “ = ^^—= ■ ( 48 ) 

dcpi Acpi dA-i AA-i 

We perform the procedure of gradient descent (40) — (48) until the integral costs for 
traffic are reduced. If further decrease in traffic is not possible, we assume, according to 
formula (34'), that new positions of N points Zj = (^t,A;), i = were found. We 

recalculate by the new values of N points the volumes of the capacities of the habitat Ui, i = 
of the corresponding Voronoi polygons. We proceed to the next step, i.e. to the 
optimization procedure (40) — (48). 

The inability to further reduce traffic means that the step of the gradient descent St 
becomes minimal and equal to S/q”, and the proportion of the variation in the traffic costs 
for one gradient descent step periodically changes sign so that the integral costs for traffic 
fluctuate near a certain constant value. 


In the procedure of gradient descent (40) — (48), in order to reduce the variability in the 
finite difference calculation of transport cost derivatives (47), (48), the relief was 
smoothed using the moving average procedure. The averaging was carried out over the 
nearest neighbours in the region of the square shape in the latitude-longitude coordinates 
with the number of points 7x7 = 49, while the original relief matrix had a resolution of 0,5° 
and was represented in the form of a matrix 360x720. 


For the initial estimate of the parameter (5 entering into the formulas (31'), (34), (34'), we 
will start with an approximate estimation of the share of transport costs in the world GDP 
in the range from 4% to 15%, depending on the development of the transport 
infrastructure. In the calculations, we will assume this share to be approximately 10%. 


We estimate the parameter (3, based on the formula (31'), assuming that the traffic costs is 
proportional to a some fraction { of the total capacity of the habitat, i.e. 


= rSi.iCl - 


^/c= 




)Ui = 


(49) 


Considering the set: the constant r, the set [u ^,..., Ujv), matrix b k = 1, —,N} and share 
we solve the equation (49) and find the coefficient J3, which should provide at the final 
stage of the process of minimizing (34'), (35') share of the traffic of the entire capacity of 
the habitat in the vicinity of 10%. 


In Figure 19,a, according to (49), a typical example of the graph of the function Tr = Tr(P') 
is constructed for N = 36,r = 0,2. The pentagram on the graph shows the value of the 
parameter (3 = 2,7292 ■ 10^, which was obtained by solving equation (49) at { = 0,11. The 
distance matrix D = [di j.i.j = was calculated from the random set of points 
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obtained according to the algorithm of Section V. Figure 19,b shows the positions of the 
random initial distribution of points (denoted by the centers of markers in the form of 
circles O) on the Earth's surface, and also the corresponding Voronoi polygons. 



0 1000 2000 3000 4000 5000 6000 


Figure 19,a. Dependence of 
transportation costs on the parameter 

P 


The Voronoi polygons for the initial distribution at N = 36 points 



Figure 19,b. The Voronoi polygons for the 
initial distribution at V = 36 points 


Figure 20,a,b shows the dynamics of the transport cost dependence on the gradient descent 
step number according to the procedure (40) — (48) with the parameters: H = 1, g = 4, 
n = 5 for two sets of points by 36 and 72 respectively. 


The dynamics of the dependence of transportation costs on the 



Figure 20,a. The dynamics of transport 
costs dependence on the step number of 
the gradient descent for N = 36 points 


The dynamics of the dependence of transportation costs on the 



Figure 20,b. The dynamics of transport 
costs dependence on the step number of 
the gradient descent for N = 72 points 


For 36 points in total, 26 steps of procedure (40) — (48) were implemented, while the cost 
of traffic decreased by 5,7% compared to the initial value. The circles in Figure 20,a,b 
indicate the transportation costs after each step of the gradient descent procedure. The 
total share of traffic costs for 36 points in relation to the entire capacity of the habitat was 
10,4%. During the whole calculation, the set of capacities of the habitat {%, ...,Uw} was 
recounted once. In Figure 20,a,b place of the recount indicated by arrows. At the final stage 
of the calculations, Figure 20,a,b shows the output of the traffic values to a certain plateau. 


Figure 21 shows the result of calculating the positions N = 36 points after applying the 
cost optimization procedure for transport (40) — (48) and one recalculation of the set of 


42 















habitats associated with each of the points. Initially, the points were chosen randomly 
according to the algorithm (21) — (24), they were marked by the center of the markers in 
the form of blue circles (O) in Figure 19,b and in Figure 21. After applying the optimization 
procedure (34'), (35'), (40) — (48) the position of the points is marked with markers in 
the form of pentagrams (★). The maxima of the displacements were 3,2° and 3,4° in 
latitude and longitude, respectively. 



-200 -150 -100 -50 0 50 100 150 200 

Figure 21. The result of optimization of transportation costs for N = 36 points, 
initially accidentally positioned on the surface of the Earth 

On the map in Figure 21 with markers in the form of crosses (X) are marked the largest 
cities^ on Earth in the number of 318. The names of 7 cities are noted, the coordinates of 
which were closest to the final positions N = 36 points. Taking into account (38) and (38'), 
the condition that the angular distance D^ 2 between a pair of points = (<^i,Ai), Z 2 = 
(<P 2 >^ 2 ) was less than one degree was chosen as a criterion for selecting cities, i.e. 0 ^ 2 < 1°, 
where z^ — coordinates of one of the cities from the set, and, and Z 2 — one of the positions 
indicated by markers in the form of pentagrams. 

The further testing of the algorithm for optimizing transportation costs (34'), (35'), 
(40) — (48) was carried out after doubling the number of points, i.e. for N = 72. Figure 
20 ,b shows the dynamics of the dependence of transport costs on the step number of the 
gradient descent. In total, 88 steps of gradient descent (40) — (48) and four recalculations 
of the set of habitats were carried out according to the formula (35'). Other parameters of 
the procedure were the following values: r = 0,2; ^ = 0,11; P = 4,7711 ■ 10^; S = 1; g = 4; 
n = 5. After the optimization procedure, the traffic costs decreased by 24,6%, the share of 
traffic relative to the entire capacity of the habitat was 8,3%. 


^The coordinates of the cities are taken from the database contained in the Matlab Mapping Toolbox. 


43 









Figure 22 shows the result of optimal positioning of Ai = 72 points. The initial positions of 
the points are denoted by markers in the form of blue circles (O), whose position centers 
are chosen randomly according to the algorithm (21) — (24). The final, optimal positions 
of the points are designated by markers in the form of pentagrams (★). The maximum 
variability of the positions of points in latitude and longitude did not exceed 12,3*^ and 
11,0°, respectively. The names of 10 cities are noted, the coordinates of which were within 
one angular degree to the total positions N = 72 points. On the map in Figure 22, the 
arrows connect the initial and final positions of points after the optimization procedure, the 
angular distance between which exceeds the value 3°. 



A 


Figure 22. The result of optimization of transportation costs for N = 72 points, 
initially accidentally positioned on the surface of the Earth 

Comparing Figure 21 and Figure 22, we can note the following. The optimal positions of the 
points (markers in the form of pentagrams) tend to fill the entire capacity of the habitat, 
that is, they tend to be distributed as uniformly as possible on the Earth's surface, taking 
into account the uneven density of the capacity of the habitat. Thus, the calculations of the 
optimal positioning of the points shown in Figure 21,22 confirm the realization of the 
minimax transport doctrine. 
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SECTION IX. 

GEOPOLITICAL CLASSIFICATION OF POINTS AND TERRITORIES 


Let us turn to the task of identifying the positioning of a set of points within the land in 
terms of marine or continental types. To do this, we define the percentage of "sea- 
continent" in global traffic, SeaContinent. Taking into account the formulas for the 
transport integrals (32), (33), we write the following relation: 

SeaContinent = 100%—-, (50) 

1+Lij=i dsea,i,]/Mj=i ‘^earth.i,; 

where dsea.ij dearth.i,; the parts of total distances, passing through the land and sea 
for a pair of points with indices i and j. 

The percentage ratio (50) is normalized as follows. Let the sums of distances passing 
through the sea be zero, i.e. dgea.i,; = 0^ then the percentage ratio is 100%, i.e. 

SeaContinent = 100%. In this case, we can talk about a purely maritime type of traffic. 
From the point of view of transportation costs, the condition that the part of the 
generalized distance be equal to zero means that transportation costs on this section are 
absent. 

In the case when the sums of distances passing over the land are zero, i.e. dearth,i,; = 
0,i,j = the percentage ratio is zero, i.e. SeaContinent = 0. In this case we are 

dealing with a purely continental type of traffic. Similar to the percentage ratio (50), 
normalized to the maritime traffic type, you can enter the percentage 

ContinentSea = 100% — SeaContinent = 100%—-“Ttwv-< (50') 

l+2]ij=l dearth,i,j72]ij=idsea,i,j 


normalized to continental type of traffic. 

The percentage ratios (50), (50') can be changed in two ways. First, by varying the 
parameters gi,g2>gz>9A> which characterize transport integrals (32), (33). Secondly, by 
varying the positions of points in relation to the traffic cost optimization procedure. The 
traffic optimization procedure was considered in the previous section. 

We rewrite the percentage ratio (50) with an explicit selection of the parameters 
gi>g 2 >g 3 >g 4 - To do this, we introduce the notation for the four transport integrals, 
4,4,4,4, which appear in (32), (33), i.e. 



Rewrite percentage (50) taking into account the transport integrals (51), then 
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1 + 5^4 S h/(di Ti h+92 S h'^93 S hy 


(52) 


where D 4 = 4 ;i,; < ^ = 1 , 2 ,3,4, 4 ;i,y is the fc-th transport integral between a pair of 

points i and j. 


Table 8 shows the results of calculating the percentage ratio (50) for the number of points 
Af = 10 and N = 10^ respectively. In ten computational experiments, the points' positions 
were generated randomly according to the procedure of Section V. 


Table 8. The results of calculating the percentage ratio of "sea - continent" for 

lArkon rt — 1 rt — t rt — 1 f»-3 n — 10-3 


SeaContinent 

#1 

#2 

#3 

#4 

#5 

SeaContinent|)v=io 

89,9% 

85,0% 

63,0% 

69,9% 

73,9% 

SeaContinent 4 = 10 ^ 

69,8% 

68,7% 

68,2% 

69,0% 

68,0% 


From the analysis of the values of the percentage ratio according to Table 8 the following 
conclusions follow. First, with the selected values of the parameters gi,g 2 > 9 z> 9 A deal 
primarily with the sea traffic type ( SeaContinent > 63%). This is not surprising, because 
the surface of the Earth, as is known, is divided into land and sea in relation to 29,2% and 
70,8% respectively. Secondly, even with the number of points exceeding N = 10^, the 
variability of the percentage ratio is sharply reduced. This means that if you do not take 
into account the specially selected positioning geometry of points, then the percentage 
ratio can be significantly changed only due to the variation of parameters g^, 92 , 93 , 94 - 

In addition to Table 8 , we consider the case when iV = 318. As a set of points, let us 
consider the database of the 318 largest cities of the Earth, mentioned in the previous 
section. After calculating the percentage (50) for them with the parameters g^ = 1 ,92 = 1. 
92 = 10“^, 94 = 10“^ it turned out to be 73,17%. The obtained value is slightly higher than 
the current estimation of the share of sea transport in global traffic equal to 2/3. Change 
the parameter ^4 so as to ensure a percentage ratio of 2/3 or 66,67%. To do this, we 
express the parameter ^4 from equation (52) through three other parameter g±,g 2 >g 3 > 
then 


^4 = ( seacon^ent “ 1)®^! + |7!^2 + (53) 

After substituting in (53) for the selected values N = 318 city points the SeaContinent = 
66,67% value as well as the quantities found using the computational experiment: 24 — 
47,9151; 24 = 30,6374; 24 = 7,6271 ■ lO'^; 24 = 5,6757 ■ lO'^, we obtain 

94 = 0,5(8,4421 ■ 10-4^1 + 5,3980 ■ 10-4^2 + 1,3438^3). (53') 

Substituting g^ = 1, 92 = 1, 93 = 10“^ into (53'), we find the modified value of the 
parameter 94 = 1,3639 ■ 10“^, which provides a percentage ratio equal to 66,67% (before 
the procedure for minimization of transportation costs). 
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By analogy with (50), we construct a percentage ratio for each point from a certain set. 
Taking into account that the matrices dsea,i,j dearth,i,;> 3re not symmetric 

in general, we will make the following percentage ratio: 


SeaContinentj = 


100 % 


l+2^l(‘^sea,i,_/+‘^sea,_/,i)/2^l(‘^earth,i,l+dearth,j,i) 


(54) 


where i = 1, The percentage ratio (54) allows us to estimate the "sea - continent" 
relationship for each point from a given set. 

Taking into account the definition in (54), we classify 318 points-cities in terms of 
preferences "sea - continent". After estimating the global percentage of 66,67% in (53), 
(53') it became clear that this value acts as a preference separator in the "sea - continent" 
relation. In other words, if for some i-th point SeaContinent; > 66,67%, then this 
circumstance is interpreted in the sense that the i-th point is oriented to the sea. 
Conversely, if some i-th point is such that SeaContinent; < 66,67%, then this means that 
the i-th point is oriented to the continent. The word "orientation" is interpreted in the 
sense of a geopolitical confrontation, that is, the "sea" strives for power over the 
"continent" and, on the contrary, the "continent" is striving to acquire the "sea". 

We note the following important circumstance in the matter of the confrontation between 
the sea and the continent. If some i-th point is focused on the sea, it means that it is rooted 
in the continent or, alternatively, is positioned continental. Conversely, if the i-th point is 
focused on the continent, it means that it is rooted in the sea, or, otherwise, has a maritime 
positioning. Thus, each point (the corresponding territory) has two attributes: 1) 
orientation and 2) positioning. From the point of view of geopolitics, both attributes with 
respect to each point have opposite names in terms of "sea - continent". This means that if 
there is a sea orientation, the positioning is continental and vice versa. The attribute 
"orientation" is characterized by a percentage ratio, SeaContinent; attribute "positioning" 
— the percentage ratio, ContinentSea. These verbal formulas should be understood in a 
metaphorical sense. All points having both a marine and continental type of orientation are 
ground-based. 

Figure 23 shows the result of a computational experiment on the calculation of percentage 
ratios (54) for each of the 318 points-cities. When calculating transport integrals, the 
following parameter values were selected: = 1, ^2 = 93 = 10“^, = 1,3639 ■ 10“^. 

Figure 23 shows a map-classification of 318 points-cities in terms of their orientation, 
either at the sea or on the continent. Points-cities are marked with markers in the form of 
pentagrams, surrounded by suitable polygons Voronoi. The cities and regions marked in 
brown are oriented to the continent. 


According to the map in Figure 23, there is a huge cluster of regions marked in blue, which 
are oriented towards the sea. This cluster is called the "Continent", it consists of Europe, 
Asia, North and Equatorial Africa, and also from a small fragment in North America. The 
continent is surrounded by a multitude of scattered territories, which consist of almost the 
whole of North America, South America, South Africa, Australia, New Zealand and Oceania. 
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Antarctica should also be included in this group. All these regions are oriented to the 
continent, i.e. they have marine positioning. In this regard, let's name the whole set of 
territories that have a sea positioning "Sea”. 


Classification of 318 points - cities in terms of orientation to the sea and the continent respectively 


9 - 0 - 



Figure 23. Classification of 318 points-cities in terms of orientation to the sea (blue) 
and to the continent (brown color) respectively 


Figures 24,a,b show individual fragments of the map in Figure 23, in which two special 
cases of mixing territories with the opposite orientation are identified. Figure 24,a shows 
North America, in which all territories are oriented to the Continent except for six ones. 
These six territories include the territories associated with cities such as Anchorage 
(69,1%), Vancouver (67,8%), Seattle (67,5%), Portland (67,1%), Edmonton (69,3%) and 
Denver (67,9%). 



Figure 24,a. Mixing the territories of the 
opposite orientation in North America 


Figure 24,b. Mixing the territories of the 
opposite orientation in England and 
North Africa 


Figure 24,b shows Europe and North Africa, in which all the territories are oriented to the 
sea except for four ones. These four territories include the territories associated with such 
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cities as Dublin (64,8%), Glasgow (64,8%), Edinburgh (64,6%) and Tripoli (66,66%). 
Percentage ratio 66,66% for Tripoli slightly less than the threshold value 66,67%, i.e. the 
area associated with Tripoli, is a border between two opposite orientations. 

Note that all of the above orientations of the territories (at sea, on the continent) should be 
interpreted from the point of view of geopolitics. The verdict about this or that orientation 
of the city (territory) was obtained by calculating the percentage (54), which was 
compared with the global percentage ratio of 66,67%. 

Let's compare the total capacity of the habitat area, oriented to the sea, Usea, with the total 
capacity of the habitat area, oriented to the continent, f/continent- calculate these values 
by the formulas: 

^Sea 2i=l,SeaContinenti>66,67% ^Continent Si=i,SeaContinenti<66,67%^i' 

where Uj, i = 1, ...,318 — the habitat capacity of each of the 318 Voronoi polygons, whose 
appearance is shown in Figure 23,24. In the sums (55) only those terms whose percentage 
ratios satisfy the corresponding inequalities were taken into account. 

After calculating the values of (55), it turned out that 100%Usea/U = 51,73%, 
100%1/continent/^^ = 48,27%, where U is the global capacity of the habitat. Thus, the 
potential for a geopolitical confrontation between territories oriented toward the sea and 
the continent is in the proportion of 51,73% to 48,27%, that is approximately equal. 

We construct analogues of clusters of territories called "Heartland” and "Rimland". These 
concepts are typical for classical geopolitics, they were introduced in the first half of the 
twentieth century by H.J. Mackinder [7] and N.J. Spykmen [12], respectively. We will also 
define the place where the notorious geographic "axis” of history passes. 

According to Figure 23 continental positioned territories (corresponding Voronoi 
polygons) are determined from the condition that the individual percentage ratio is more 
of 66,67%. Let such territories be N^. Let us find among these territories the average value 
of the percentage ratio C according to the following formula: 

c = ^ Sfii,seaContinenti>66, 670/0 SeaContinent j. (56) 

The presence of the average percentage ratio calculated according to (56) allows us to 
divide the territories into two categories, in which the individual percentages are in two 
ranges: [66,67%; C^ and \C; max SeaContinent,]. Territories falling within the first range 

l<i<318 

we refer to Rimland, and within the second range — to Heartland. 

As candidates for the geographical location of the "axis” of history, using the weighted 
average, we define two points. The weights will be the corresponding values of the 
capacities of the habitat of each of the territories. 
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The first point acts as a center of "gravity" for all territories oriented to the sea or 
continental-positioned. The position ((pc,^c) of this point is found according to the 
formulas: 


V318 mil- ySlS t 

_ '^i=l,SeaContinentj>66,67 %^ _ '^i=l,SeaContinentj>66,67%^i^i 

~ y318 ~ ^ ~ y318 ~ 

-^i=i,SeaContinentj>66,67% -^i=i,SeaContinentj>66,67% 


(57) 


The second point acts as the center of gravity for all territories attributed to Heartland. The 
position {(pfj, Xfj) of this point is found by formulas: 


_ 2]i=i,seaContinent(>C ^ _ 2]i=i,seaContinentj>C'^i^i 

Vh ~ y318 ~ > ~ y318 ~ 

-^i=l,SeaContinentj>C -^i=l,SeaContinentj>C 


(58) 


In the framework of the above model considerations, we will construct a map with 
Rimland, Heartland and a pair of points (57), (58), applicants for the geographical "axis” of 
history. On Figure 25 the result is shown. Black colour is allocated to the territories 
attributed to Heartland, red colour is allocated to the territories attributed to Rimland. 



Figure 25. Positioning of Heartland, Rimland and a couple of geographic “axes” of 

history 


It turned out that the division of territories into Rimland and Heartland is determined with 
the value of the percentage ratio equal to 73,7%, because C = 12,1%. The Heartland has hit 
83 of the territory, and in the Rimland — 119, total Nc=202. In terms of the capacity of the 
habitat, the shares of Heartland and Rimland in the Continent were 47,52% and 52,48% 
respectively. 
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Finally, two points representing the geographical "axis” of history were calculated and 
plotted on the map (two markers in the form of rhombuses). Larger rhombus indicates the 
point with coordinates (pc = 31,8386°, Ac = 42,9602°, which correspond to the city of 
Karbala in Iraq. Smaller rhombus denotes the coordinates (p^ = 39,1497°, = 87,7375°, 

which correspond to the nameless positions in the Xinjiang Uygur Autonomous region of 
China. 

Due to the importance of the terms "Heartland” and "Rimland” from the point of view of 
geopolitics, we will define within Heartland and Rimland a subdivision into two groups of 
regions relative to the corresponding average values of percentage ratios. In other words, 
we divide Heartland into Heartland 1 and Heartland 2 and, similarly, Rimland to Rimland 1 
and Rimland 2. 


The Heartland division in Heartland 1 and Heartland 2 



Figure 26. The Heartland division in Heartland 1 and Heartland 2 


Figure 26 shows the map on which the desired subdivision into two Heartland is built, with 
Heartland 1 having a higher percentage ratio SeaContinent than Heartland 2. On the map, 
the boundaries of the territories classified as Heartland 1 are painted black, the boundaries 
of the territories classified as Heartland 2 are painted dark red. 

The study of the map in Figure 26 leads to unexpected conclusions for classical geopolitics. 
Inside the Continent there is a compact core, called Heartland 1, which is surrounded by a 
belt called the Heartland 2. To Heartland 1 we can include Central Asia, Tibet, most of 
China, part of the Russian Federation. 

Figure 27 shows the map on which the desired subdivision into two Rimlands is built, with 
Rimland 1 having a higher percentage ratio SeaContinent than Rimland 2. On the map, the 
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boundaries of the territories classified as Rimland 1 are painted red, the boundaries of the 
territories designated as Rimland 2 are colored magenta. 


Subdivision Rimland on Rimland 1 and Rimland 2 



Figure 27. Subdivision Rimland on Rimland 1 and Rimland 2 


Positioning of territories classified as Rimland 3 and Rimland 4 



Figure 28. Positioning of territories classified as Rimland 3 and Rimland 4 


The study of the map in Figure 27 shows a strong non-connectivity, both Rimland 1 and 
Rimland 2. In Rimland 1 reveals two large, compact clusters, is Central Europe and Central 
Africa. 
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It remains to examine in more detail the internal structure of territories oriented to the 
continent, i.e. positioned as marine ones. For these territories, the individual percentage 
ratio is less than or equal to 66,67%. We divide these territories into two classes with 
respect to the corresponding average value of the percentage ratios. Let's name these two 
classes of territories "Rimland 3" and "Rimland 4”. At the same time, the percentage ratios 
of the Rimland 3 territories are larger than the percentages of Rimland 4. 

Figure 28 shows the map on which the desired division was built into two classes: 
Rimland 3 and Rimland 4. On the map, the borders of the territories classified as Rimland 3 
are painted red, the boundaries of the territories referred to the Rimland 4 category are 
painted in magenta colour. 

Attention is drawn to the map in Figure 28 attributing the USA, most of Canada, the middle 
part of South America, and also South Africa to Rimland 3. Finally, Central America, the 
south of South America, Australia and most of the Antarctic are parts of Rimland 4. 



Figure 29,a. Positioning of a point and 
territory that has the maximum 
percentage ratio 



•aa ■ * - i I . i - 
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Figure 29,h. Positioning of a point and 
territory that has a minimum percentage 
value 


We will find a couple of territories from our list that correspond to the maximum and 
minimum values of the percentage ratio. Figure 29 shows the result. Figure 29,a shows the 
city of Lanzhou and the corresponding territory, which have a maximum value of a 
percentage ratio equal to 82,84%. Figure 29,b shows the city of Auckland and the 
corresponding territory, which have a minimum percentage of 39,15%. The given cities and 
territories from the selected list are maximally removed from each other in the geopolitical 
sense of the word, with the city of Lanzhou located in Heartland 1 and the city of Oakland in 
Rimland 4. 


Table 9. The final geopolitical classification of points and territories 


Nomenclature 
of classes 

Sea 

Continent 

Rimland 

Heartland 

Rimland 4 

Rimland 3 

Rimland 2 

Rimland 1 

Heartland 2 

Heartland 1 

SeaContinent, % 

139,15:56,651 

[56,65; 66,671 

[66,67:70,811 

[70,81:73,741 

[73,74:77,961 

[77,96:82,841 

U = 7,64 

1,5865 

2,1010 

0,9404 

1,1339 

0,9538 

0,9244 


Put together our classification of points and territories in terms of geopolitics. The result is 
presented in the form of a Table 9. 
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In Table 9 there is collected the nomenclature of 6 classes of points and territories, named 
in terms of geopolitics. The intervals of the values of the percentage ratio of points 
(territories) that characterize a particular class are given. The last line of Table 9 shows the 
capacities of the habitat of each class, as well as the sum of all individual capacities, equal to 
the global capacity of the habitat U. 

Taking into account Table 9, we will build a map on which we will draw the geopolitical 
lines that separated Heartland 1, Heartland 2, Rimland 1, Rimland 2, Rimland 3, Rimland 4. 
In other words, we will construct a map representing the geopolitical structure of the 
Earth's surface according to this model. The result is shown in Figure 30. 



A 

Figure 30. The geopolitical structure of the earth's surface 


Compatible map of the geopolitical structure of the Earth, shown in Figure 30, with the 
current political map of the world. This will make it possible to build a list of States through 
whose territory certain geopolitical lines pass. These lines can be interpreted as "fracture” 
taking into account all the current and, perhaps, the past history of an individual state of 
the list. Since it is extremely difficult to combine all these lines on one map, we will 
consider some fragments of the global map. 

Build a separate Heartland plus associated countries borders. Figure 31 shows the result. 
After studying this map, you can conclude the following. Russia is almost completely 
included in Heartland, more precisely in Heartland 2, there are only small intersections 
with Rimland. China is almost completely included in Hartland 1. 
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Heartland plus the state borders of the associated countries 



40 60 80 100 120 140 160 180 
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Figure 31. Heartland plus the state borders of the associated countries 

A line dividing Rimland 1 and Heartland 2 passes through Ukraine from north to south. 
Taking into account current and past events in Ukraine, it can be concluded that the line 
dividing Rimland 1 and Heartland 2 is the line of the geopolitical fracture. 


Political map of Europe plus geopolitical marking 



Figure 32. Political map of Europe plus geopolitical marking 

It is interestingly that this same line crosses Turkey, and in Rimland 1 there is somewhere 
about a third of the territory of Turkey, and in Heartland 2 — two-thirds. This line passes 
through a part of Syria. It turns out that Oman belongs to Heartland 2. The same fracture 
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line passes through the territory of Thailand, cutting it almost in half. It turns out that 
Taiwan refers to Rimland 1. With the Korean peninsula, it is especially interestingly, two 
lines pass along it: the first line separates Heartland 1 from Heartland 2, the second line 
separates Heartland 2 from Rimland 1. Note that the first line is pretty close to the line 
dividing the two Koreas. Finally, it is rather interestingly, Japan almost in half to share the 
line separating Rimland 1 and Rimland 2. This line is not considered as a geopolitical 
fracture line. 

Compatible political map of Europe and geopolitical lines defined on Figure 30. Figure 32 
shows the desired map. 

The visual inspection of the map in Figure 32 shows that within Europe there is a nucleus, 
called Rimland 1, surrounded by territories called Rimland 2. The line separating 
Rimland 1 and Rimland 2 most significantly crosses the territories of France, Germany, 
Poland and Finland. Finally, England is almost equally divided between Rimland 2 and 
Rimland 3. 

In the terminology adopted above, we will continue to study the issue of the presence 
(absence) geopolitical fracture lines. Let's pay attention to Africa. Let us gather together the 
political map of Africa and draw geopolitical line, a list of which is shown in Figure 30. 
Figure 33 shows the final map for Africa. 



Figure 33. Political map of Africa plus geopolitical marking 


The general examination of the map in Figure 33 shows that the State borders of individual 
countries on the African continent and geopolitical lines are not significantly coordinated. A 
more detailed study of the map shows the following. Somalia is divided in half by the lines 
separating Rimland 1 and Rimland 2. Given the past and current within Somali 
circumstances, we can consider the line separating Rimland 1 and Rimland 2 as a 
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geopolitical fracture. The Democratic Republic of the Congo, Algeria, the Libyan Arab 
Jamahiriya unite territories from three geopolitical classes: Rimland 1, Rimland 2 and 
Rimland 3. Through Gabon, Congo, Cote d'Ivoire, Guinea, Mauritania, Egypt, is the line 
separates Rimland 1 and Rimland 2, and through Angola — Rimland 2 and Rimland 3. 

In connection with the analysis of the location of geopolitical lines in the North and South 
America [Figure 30], we note that there is no reason to consider them geopolitical fracture 
yet. If something negative in the geopolitical sense of the word will happen, so it is likely to 
be relevant to the line dividing Rimland 2 and Rimland 3 in northwest North America. The 
other geopolitical lines are not fracture, including because the neighbouring territories 
have a single marine positioning. 
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SECTION X. 

THE RESULTS OF THE OPTIMIZATION OF TRANSPORT COSTS 


Let us pass to the final task of the transport costs minimizing for the number of points 
equal to N = 318. As initial values for the 318 points, we choose the positions of the largest 
cities from the available database. Let us study the geometry of the optimal positioning of 
points within the land in the context of the "sea - continent” dichotomy. In this connection, 
we solve three problems of optimal positioning of points: 

1 ) purely marine type of transportation of goods, when the cost of goods transportation by 
sea is zero, SeaContinent = 100%; 

2 ) purely continental type of transportation of goods, when the cost of transportation of 
goods by land is zero, SeaContinent = 0%; 

3) mixed type of transportation of goods, when the cost of transporting goods by sea and 
by land is comparable, for example, following (52) — (53'), we assume that 
SeaContinent = 66,67%. 

Note that, from the point of view of common sense, the first and second tasks are 
impossible in reality, since the transportation costs are minimized, but always remain non¬ 
zero for any type of transport. However, the solution of these problems is important for 
understanding the solution of the third problem in terms of the "sea - continent” 
dichotomy. According to (52), the condition for zeroing the transportation costs by sea is 
reduced to the condition ^4 = 0 , similarly, zeroing the transportation costs by land is 
reduced to the conditions: g\ = § 2 = 93 = 0 - 

In order to present the solutions of the first and second tasks in general terms, it is 
necessary to take into account the minimax transport doctrine introduced earlier. Let the 
transport costs by the sea equal zero, in this case it is necessary to present the geographical 
map of the world without oceans, that is, only the land remains, within which, according to 
the minimax transport doctrine, the points are maximally removed from each other, taking 
into account the heterogeneity of the habitat density distribution. If transportation costs by 
land are equal to zero, it is necessary to present a geographical map of the world without 
continents, then for positioning of points there is only a coastline over which they are 
evenly distributed taking into account the heterogeneity of the habitat capacity density. 
Thus, the maritime type of transport communications promotes the movement of points 
into the continents and, conversely, the continental type of transport communications 
moves the points closer to the shoreline. The presented conclusions in a brief metaphorical 
form can be reduced to the formula of interpenetration of the sea and land: "sea ^ 
continent”. 

Let us solve the first problem, when it is considered that the goods transportation costs by 
sea are zero, i.e. SeaContinent = 100%. Figure 34 shows the total optimization of 318 
points in the first task. As the parameters for the optimization procedure discussed in 
Section IX, the following values were chosen: r = 0,2; ^ = 0,11; P = 1,7753 ■ lO'’^; S = 1; 
q = 4; n = 5, g± = 1, 52 = 1< 93 = 10“^, g^ = 0. Initially, the positions of the points 
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coincide with the positions of 318 cities. To achieve the optimum within the framework of 
this procedure, it took 30 steps and one recount of the habitat capacities for each of the 
points. 


The total cost optimization for transportation for A^ = 318 points with 
a percentage ratio of SeaContinent = 100% 



Figure 34. The total cost optimization for transportation for iV = 318 points with 
a percentage ratio of SeaContinent = 100% 


In Figure 34, the centers of the markers in the form of blue circles (O) are the initial 
positions of the points. The centers of markers in the form of pentagrams (★) indicate the 
positions of points after transportation costs minimizing. As a result of the optimization 
procedure, the total transportation costs decreased by 39,6% and amounted to 6,6% of the 
total capacity of the habitat. In Figure 34, the names of cities marked the pentagrams in the 
number of 13, whose positions were within 0,5° angular distance from one of the cities of 
the original list. The arrows indicate the transitions from the initial positions to the optimal 
ones when the angular distance exceeds the value 1,5°. 

The study of the optimum positions of the points on Figure 34 demonstrates 
unambiguously the movement of points into the continents, which corresponds to the 
action of the minimax transport doctrine. 

Let us proceed to the solution of the second problem, when it is assumed that the costs of 
transporting goods over land are zero, i.e. SeaContinent = 0%. Figure 35 shows the total 
optimization of 318 points in the second task. The following values were chosen as 
parameters of the optimization procedure: r = 0,2; ^ = 0,185; P = 2,5554 ■ lO"*^; E = 1; q = 
4; n = S, = 0, g 2 = 0, gs = 0, ^4 = 10“^. Initially, the positions of the points coincided 
with the positions of 318 cities. To achieve the optimum within this procedure, it took 40 
steps and three recalculations of the habitat capacities of each of the points. 
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In Figure 35, the centers of the markers in the form of blue circles (O) are the initial 
positions of the points. Centers of markers in the form of pentagrams (★) indicate the 
positions of points after the transportation costs minimizing. As a result of the optimization 
procedure, the total transportation costs decreased by 9,4% and amounted to 16,8% of the 
total capacity of the habitat. In Figure 35, the names of cities marked pentagrams in the 
number of 57, whose positions were within 0,5° angular distance from one of the cities of 
the original list. 


The total cost optimization for transportation for A = 318 points with 



Figure 35. The total cost optimization for transportation for iV = 318 points with 

a percentage ratio of SeaContinent = 0% 

The study of the optimum positions of the points in Figure 35 demonstrates the obvious 
movement of points closer to the shoreline in accordance with the action of the minimax 
transport doctrine for this case. Unfortunately, the available computational resources did 
not allow to use the optimization algorithm of the previous section with more advanced 
values of the parameters (iV, E, q, n), which would ensure the movement of all points to the 
shoreline. 

In the context of the general description of the solutions of the first and second tasks, the 
meaning of the solution of the third problem, which has a direct relation to reality, becomes 
clear. In the case of non-zero transportation costs, both by sea and by land, the distribution 
of points is similar to something in between maritime and continental types of the 
positioning points of the first two tasks. As an example, consider the intermediate type, 
when the percentage ratio is chosen equal to SeaContinent = 66,67%, and the transport 
parameters according to (52) — (53') — .^i = 1. 52 = 1< ds = 10“^, 54 = 1,3639 ■ 10“^. 

Figure 36 shows the total transportation cost optimization for N = 318 points, the initial 
location of which coincided with the position of the largest cities from the existing 
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database. As parameters of the optimization procedure discussed in Section IX, the 
following values were chosen: r = 0,2; ^ = 0,11; P = 8,2564 ■ 10^; S = 1; g = 4; n = 5. To 
achieve the optimum within the framework of this procedure, 27 steps and two 
recalculations of the habitat capacities of each of the points were required. In Figure 36, the 
centers of the markers in the form of blue circles (O) are the initial positions of the points. 
Centers of markers in the form of pentagrams (★) indicate the positions of points after the 
transportation costs minimizing. As a result of the optimization procedure, the total 
transportation costs decreased by 7,8% and amounted to 10,1% of the total capacity of the 
habitat. In Figure 36, the names of cities marked pentagrams in the number of 31, whose 
positions were within 0,5° angular distance from one of the cities of the original list. 

The total cost optimization for transportation for 318 points with 



Figure 36. The total cost optimization for transportation for iV = 318 points with 
a percentage ratio of SeaContinent = 66,67% 

As a result of solving the three problems of optimal transport positioning of a set of points 
mentioned above, we can formulate the following general problem. To choose the transport 
cost functional, so that the optimal positions of some set of points in relation to this 
functional correspond to the real positions of the cities acting as logistic centers. The 
computer resources available to the author at the current time are scant to solve this 
problem. 
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CONCLUSION 


In this article we construct the mathematical model of geopolitics. The basis of the model is 
the central concept, called the "capacity of the habitat”. This concept is defined and 
calculated. The function of density of habitat capacity is constructed, which depends on the 
distribution of mean annual temperature and precipitation on the Earth's surface. A high 
correlation between the density of the capacity of the habitat and the population density of 
the Earth was established. 

The density of habitat capacity has been calculated for each State. Among the record 
holders in the diminishing sequence were expected: Russia, the USA, Brazil, China, 
Australia, etc. An indicator has been determined and studied that has a specific habitat 
capacity per capita. The ranking of countries and territories by this indicator has been 
made. Particular attention is paid to the ratio of these indicators of individual countries in 
relation to the Russian Federation. 

Studied the issue of the relationship of the density capacity of the habitat and topography 
was studied. The areas where 50% of the population and 50% of the habitat are 
concentrated are built and compared. The fields of the density gradient of the habitat are 
being constructed. 

Classification of countries and territories is carried out in terms of "high - low" and 
"favorable - unfavorable”, i.e. in four categories, taking into account the relief and the 
density of the habitat. Maps of territories of all four types are constructed. An index of the 
diversity of individual territories and States is introduced and calculated. 

An algorithm for random drawing of points on the surface of the land is developed, taking 
into account the density of the capacity of the habitat. This algorithm is necessary to 
develop a solution to the optimal positioning problem of points in terms of minimizing the 
transport costs between them. 

A mathematical model of transport flows between individual points, which act as logistic 
nodes, is formulated. As the basis were used the well-known gravitational model and the 
exponential type of dependence on generalized distances between points. A minimax 
transport doctrine is introduced, which ensures, when optimizing the positions of the 
points, a minimum of transportation costs with a maximum of occupancy of the habitat 
taking into account the heterogeneous distribution of the density of the habitat capacity. 

The method of calculating the generalized distances between points is described in detail, 
taking into account the difference in the sections passing over land and sea. These 
distances are estimated in terms of the specific energy costs for moving the weight unit of 
the conditioned cargo. 

An algorithm for minimizing transport flows for an arbitrary set of points is constructed 
and tested. This algorithm is some modification of the gradient descent procedure taking 
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into account the shoreline and the non-connectivity of the domain of determining the 
positions of the points of the transport cost functional. 

Within the calculation of global traffic, a special indicator has been constructed, called the 
percentage ratio "sea - continent”. Based on this indicator, the points (territories] are 
classified in geopolitical terms. This indicator allowed us to formalize the introduction of 
such well known in geopolitics concepts as Heartland and Rimland. 

The division of Heartland into Heartland 1 and Heartland 2 was implemented. It turned out 
that the core of Hartland called Heartland 1 consists of Central Asia, Tibet, most of China 
and part of the Russian Federation. If Heartland and Rimland in geopolitics are usually 
referred to the "Continent", then everything else is attributed to the "Sea”. It turned out that 
the capacities of the habitat of the territories positioned in the Sea and the Continent is 
about the same. The positioned in the Sea territories were divided into two groups with 
respect to the average value and named Rimland 3 and Rimland 4. The cities (Auckland and 
Lanzhou] are identified from the list under consideration, consisting of 318 cities, with 
minimum and maximum percentage ratios. A complete classification of all the city points 
and the corresponding areas from the adopted list in terms of geopolitics is carried out. 

Combined global and regional maps containing political and geopolitical markings are 
constructed. The analysis of these maps for the presence of so-called geopolitical "fracture” 
is made. The presence of fracture can be said in the event that some geopolitical lines pass 
not along the border of an individual State, but through its territory, penetrating it deeply 
into it. 

Within the framework of solving the minimizing transportation costs problem, 318 points 
were considered, the initial location of which coincided with the position of 318 largest 
cities on Earth. In a single complex, three subtasks were solved numerically. When solving 
each of the three subtasks, the minimax transport doctrine was essential. In the first 
subtask, it was assumed that the costs of transport by sea were zero; in the second task it 
was considered that the costs for transport by land were zero and, finally, in the third 
subtask, the real costs of transportation were chosen. It is shown that the optimal 
configuration of the points of the third subtask represents something in relation to the 
optimal configurations of the first two subtasks. 
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Mathematical Model of Geopolitics 


The mathematical model of geopolitics is a conditional name for several models, which 
are naturally connected and act as accompaniment to the main theme — geopolitics. All 
constructed work models made the transition to computing experiment, the results of 
which are given and discussed. The central concept of the mathematical model of 
geopolitics is introduced - the capacity of the habitat. Geopolitics is the climate, the 
relief, the logistics features of global commodity flows, the geopolitical confrontation in 
terms of the “sea-continent”, i.e. all that constitutes the material complex of living 
conditions of the inhabitants of the Earth. This complex, to a large extent, mediates the 
population's behavior from the political point of view. The author does not adhere to the 
position of natural determinism, which acts in the form of geopolitics, but tries to 
delineate the scope of the manifestation of geopolitics in real politics. The paper provides 
clarification and generalization of the generally accepted geopolitical classification of 
territories in terms of orientation and positioning either at sea or on the continent. The 
mathematical model of transport costs minimization for an arbitrary number of points 
acting in the form of logistic centers is formulated. 
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